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DENSITY OF LIQUID HELIUM-3 BETWEEN 0.045 AND 1.3°K* 


John E. Rives and Horst Meyert 
Department of Physics, Duke University, Durham, North Carolina 
(Received August 7, 1961) 


We have measured the density of liquid helium- 
3 at various pressures up to 28.0 atm in the tem- 
perature range between 0.045 and 1.3°K. At pres- 
sures above about 15 atm, an unexpected minimum 
in density was found to occur, the temperature of 
this minimum increasing to 0.09°K when the pres- 
sure is increased to 28 atmospheres. For a pres- 
sure of 0.18 atm and below 0.08°K the coefficient 
of expansion was found to be proportional to T, as 
predicted by Brueckner.’ 

The method used for these measurements was 
the same as that used by Lee, Fairbank, and 
Walker,” who determined the density of liquid 
helium-3 down to about 0.15°K. It consisted in 
measuring the dielectric constant of the liquid 
from the resonant frequency of an LC circuit, in 
which the capacity was a coaxial cylindrical ca- 
pacitor containing the liquid. The density of the 
liquid was then obtained using the Clausius-Mo- 
sotti relation. The assumption was, as in earlier 
work,? that the polarizability stayed constant with 
temperature and pressure. 

The cylindrical capacitor was made of parts 
machined from electrolytic copper, the gap being 
4x10~? mm between the concentric cylinders. 
These were held together at their upper extremity 
by an araldite seal capable of withstanding pres- 
sures at least up to 29 atmospheres. Helium-3 
could be introduced into the gap through a 2 x107?- 
mm i.d. capillary. The inductance of the tank 
circuit consisted of a coil of copper wire soldered 
directly on top of the capacitor and shielded from 
external pickup by a thin copper cylinder connect- 
ed to the ground side of the capacitor. The tank 
circuit was thermally connected via 500 No. 34 











































copper wires to the paramagnetic salt pill made 
up of about 30 g of freshly grown chromium 
methylammonium sulphate mixed with Apiezon 

oil J. The salt was adiabatically demagnetized 
from a field of about 12 kilogauss at about 0.45°K, 
and the lowest temperatures reached were about 
0.035°K. The temperature was measured by two 
470-ohm Speer carbon resistors, the first sol- 
dered to the tank circuit, the other embedded in 
the salt pill. These resistors were calibrated in 
each series of experiments against the suscepti- 
bility of the chrome alum, which in turn had been 
directly measured against that of cerium magnesi- 
um nitrate down to 0.06°K in a different apparatus. 
The temperatures below 0.06°K were determined 
by extrapolation of the calibration curve of the car- 
bon resistors. While the temperature should be 
known to within 6% down to 0.07°K, the systematic 
error may increase to about 15% at the lowest 
temperatures. 

The tank circuit was part of a Franklin oscil- 
lator® and the frequency was measured by a Hew- 
lett Packard 524B Electronic Counter. The reso- 
nance frequency of the tank circuit was about 31 
Mc/sec and was stable within 5 cycles over peri- 
ods of an hour. The power input into the tank cir- 
cuit was so small that it did not affect the warming 
up of the sample even at the lowest temperatures. 
The change in frequency between 0.045°K and 1.2°K 
due to the change in density of liquid helium-3 
at constant pressure was only of the order of 4 
kc/sec, and hence during the experiments, the 
pressure had to be held very constant, as a change 
in pressure of 0.1% produced a frequency shift in 
the tank circuit of the order of 30 cps. During 
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FIG, 1. (a) Density of helium-3 normalized to its value at 1.2°K for several pressures. --- Data of Lee, Fair- 


bank, and Walker. (b) Some density curves of liquid helium-3 below 0. 15°K. 


each series of experiments, the frequency of the 
empty tank circuit was first measured over the 
whole temperature range. At several tempera- 
tures above 1°K the density of liquid helium-3, 
obtained in arbitrary units, was normalized to 
the accurate density values obtained by Sherman 
and Edeskuty* and the results were found to be 
very consistent. 

Our results are shown in Fig. 1(a), where for 
convenience we have plotted the ratio p(T, P)/ 
p(1.2,P), where p(1.2, P) is the density at 1.200°K 
and at the pressure P as found by Sherman and 
Edeskuty.* Our density results, extrapolated to 
zero pressure, agree satisfactorily with those of 
Kerr and Taylor® and are somewhat different from 
those of Lee, Fairbank, and Walker.” As can be 
seen from Fig. 1(a) and especially Fig. 1(b), the 
density at pressures above about 15 atm and at 
temperatures below 0.1°K passes through a mini- 
mum, the temperature of this minimum increas- 
ing with pressure. Several experiments were 
carried out to check this effect and gave re- 
producible results. The coefficient of expansion 
ap = (1/V)(8V/8T)p as derived from the p values 
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is plotted in Fig. 2 for a few pressures. In the 
region below 0.08°K, we found ag 1g atm 

= -(0.12+0.02)T in qualitative agreement with 
Brueckner’s theory’ which predicts Ap = -0.0767. 
The figure also suggests that the density minimum 
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FIG. 2. The expansion coefficient of liquid helium-3 
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may be found for pressures below 15 atm at tem- 
peratures below 0.045°K, As a consequence, be- 
low 0.09°K, the entropy of isothermal compression 
AS, where 


P 
AS=S_-S = [ Va,dp, 
P  svp svp bd 


will pass through a maximum. Here Ssvp is the 
entropy at saturated vapor pressure and V is the 
molar volume. 

The quantitative value of AS as obtained from 
our density measurements can be checked by 
specific heat measurements at various pressures. 
At temperatures above 0.3°K, the agreement with 
AS as calculated by Brewer and Daunt® is very 
satisfactory. On the other hand, there are some 
systematic differences between our curves and 
those calculated from the earlier density meas- 
urements of Lee, Fairbank, and Walker.” It may 
be mentioned that no discontinuity of a», charac- 
teristic of a second order transition to another 
phase, such as predicted by Emery and Sessler’ 
and by Brueckner et al.® for liquid helium-3, was 
found in these experiments over the whole tem- 
perature and pressure range. Previous specific 
heat measurements® on liquid helium-3 under a 
pressure of 0.18 atm down to 0.008°K had already 
failed to show this transition. 

The temperature at which the maximum of den- 
sity occurred for a given pressure agreed very 
well with the results of Brewer and Daunt® up to 
about 15 atmospheres and at higher pressures 
agreed better with those of Sherman and Edeskuty.* 





A new and improved apparatus is being con- 
structed in order to extend our measurements 
down to 0.01°K and obtain more accurate values 
for a, below 0.07°K. 

The initial stages of this experiment were super- 
vised by Dr. W. M. Fairbank, whom we wish to 
thank for his valuable advice. We also acknowl- 
edge the help of Mr. C. Boghosian in some exper- 
iments. We also thank Dr. N. Kurti, Dr. H. Lon- 
don, and Dr. D. M. Lee for helpful discussions. 


*This work has been supported by the National Science 
Foundation, the Alfred P. Sloan Foundation and the 
Office of Ordnance Research. 
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SELF-DIFFUSION COEFFICIENT AND NUCLEAR SUSCEPTIBILITY OF LIQUID HELIUM-THREE* 


A. C. Anderson, W. Reese, R. J. Sarwinski, and J. C. WheatleyT 
Department of Physics, University of Illinois, Urbana, Illinois 
(Received August 14, 1961) 


Recent measurements of the specific heat’’? and 
the thermal conductivity® of liquid He® at low pres- 
sures have demonstrated, for temperatures less 
than approximately 0.04°K, the temperature de- 
pendences predicted by the Fermi liquid theory 
of Landau* and Abrikosov and Khalatnikov.* How- 
ever, the magnetization self-diffusion coefficient 
measured using the spin-echo method by Ander- 
son, Hart, and Wheatley® (AHW) obeyed a T~% 
temperature dependence down to 0.03°K in con- 
tradiction with the T~* dependence predicted by 
Hone’ on the basis of the Fermi liquid theory. 

As pointed out by AHW, their experiment was 
susceptible to error, particularly at the lowest 
temperatures, because of long equilibrium times 
and difficulties with thermometry. Consequently, 
the experiment was repeated under much improved 
conditions and extended to lower temperatures. 
We now find the diffusion coefficient to be propor - 
tional to T~ below 0.04°K -0.05°K. The nuclear 
susceptibility is constant below 0.1°K and equal 

to the value previously measured (AHW) within 
experimental accuracy. The lowest temperature 
of these measurements was 0.02°K. The meas- 
urements were carried out at a pressure of 12 

cm Hg. 

The cell containing the He* was made from Epi- 
bond 100A,° and the plug which closed the cell 
was sealed with Epibond 121.* Both resins are 
nonmagnetic.® The space containing He*® consisted 
of a cylinder 1.2 cm high and 1.2 cm in diameter 
with a coaxial cylindrical tail 5.1 cm long and 
0.25 cm in diameter. The inside surface of the 
cell was lined with 400, 0.004-inch diameter, 
99.999% pure annealed copper wires about half 
of which did not extend into the tail. The larger 
cylinder contained 1.29 g of powdered cerium mag- 
nesium nitrate (CMN) which served as a thermom- 
eter. The CMN was prevented from entering the 
tail by a cotton plug. The He® was introduced into 
the cell by a 1/64-in. 0.d. by 0.003 -in. wall 70-30 
cupro-nickel tube which had a negligible magnetic 
effect.° A pi-wound, 3000-turn pickup coil of 
0.002-in. diameter manganin wire was placed 
over 2.0 cm of the tail at the greatest possible 
distance from the CMN. The above construction 
gave internal equilibrium times of 15 minutes or 
less at the lowest temperatures while providing 
sufficient separation between the He* being meas- 
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ured and the CMN to give a negligible temperature. 


dependent magnetic gradient. The He* was cooled 
by 37.8 g of potassium chrome alum using tech- 
niques described elsewhere.*° 

The thermometer was calibrated in the temper- 
ature range of 1°K to 0.3°K using a 470-ohm Speer 
resistor which previously had been calibrated in 
this temperature range with a CMN magnetic ther- 
mometer. The accuracy of the temperature scale 
is estimated to be 2%. 

The He® used had a He* content of less than 3 
parts in 10° and a tritium decay rate of less than 
900 disintegrations/sec. The apparatus outside 
the cryostat and the method used to obtain the ech- 
oes were similar to those described elsewhere.” 

The self-diffusion coefficient, D, was meas- 
ured by the method of spin echoes using an rf 
field of frequency 85 kc/sec and of amplitude 
(2H,) 0.6 gauss measured using a nutation tech- 
nique. A 90°-180°-180° pulse sequence was used 
with a fixed time delay between the 90° pulse and 
the first 180° pulse and a variable time, 7, be- 
tween the two echoes. The resulting echoes were 
displayed on an oscilloscope and photographed. 
The ratio, R, of the amplitudes of the two echoes 
is given by R = exp(-7G"Dr*/12 -7/T,), where G 
is the magnetic field gradient, y = 2.038 x 10* 
(gauss-sec)™, and T, is the transverse relaxa- 
tion time. The self-diffusion coefficient was 
calculated neglecting the effect of T,, which had 
been found to be greater than 0.74 sec near T 
=0.03°K.° Measurements of D were made with G 
alternately parallel and antiparallel to the steady 
field so as to eliminate to first order the effects 
of residual gradients. 

The susceptibility was measured with the same 
apparatus but with no applied gradient and with 
only one 180° pulse applied between 2.5 and 5 
msec after the 90° pulse. No systematic effect 
of the time of the echo on the echo height was 
observed, even at the lowest temperatures. In 
order to increase the echo height, the He* was 
magnetized in a field of about 100 gauss which 
was turned off just before the pulse train was 
initiated. This field remained the same through- 
out the experiment. The gain of the amplifiers 
was measured at each point and the echo heights 
were corrected to constant gain. With the above 
conditions, the echo height, corrected for gain, 
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FIG, 1. 


Self-diffusion coefficient, D, vs 1/T. The 


points represent averaged data. The curve has a slope 


of 2.00 +0.02 for T less than 0.04°K. The curve la- 
beled AHW is the smoothed curve from reference 6. 


FIG, 2. 


Effective magnetic 


temperature T*=(x7)1, 153°K/X 


vs 1/T, 


(xT)1,153°K is the 


product of the susceptibility and 
temperature at 1.153°K. x is 
the susceptibility at temperature 
T. The points within each tem- 
perature interval represent 


averages for each run, 
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should be proportional to the nuclear susceptibil- 
ity. 

The results of the diffusion measurements are 
shown in Fig. 1. For T less than 0.04°K to 0.05°K 
we find that D =1.54x10-°T~ cm? °K?/sec. The 
exponent of T in this temperature range is 2 with 
an estimated accuracy of 1%. The coefficient of 
D has an estimated accuracy of 7% on the basis 
of uncertainties in the temperature, time scale, 
and gradient. At higher temperatures the results 
are almost parallel to those of AHW but displaced 
about 5% toward lower T. Their low-temperature 
points lie below ours, giving a slower temperature 
dependence than that observed by us. The temper- 
atures used by AHW were made uncertain, espe- 
cially at the low-temperature end, by the presence 
of spurious magnetic materials. Moreover, their 
equilibrium times were several hours at the low- 
est temperatures, indicating that very small spu- 
rious heat inputs could adversely affect the re- 
sults. 

In Fig. 2 we have plotted (x7)1.153°K/x vs 1/T, 
where the numerator is the product of the echo 
height at the bath temperature, the bath tempera- 
ture (1.153°K), and the ratio of the density of the 
liquid at 0°K to that at 1.153°K as estimated from 
the data of Lee, Fairbank, and Walker’; and 
where the denominator is the echo height at the 
value of 1/T given on the abscissa. If it is as- 
sumed that He* obeys Curie’s law at the bath 
temperature, then the ratio (x7)1,153°K/x is the 
effective magnetic temperature, T*, of the He® 
at low temperatures. Below 0.1°K we find T* to 
be constant within experimental accuracy. The 
average of the data for run 4, which is the best 
run for this purpose since the check at the bath 
temperature was made at its conclusion, gives 
T* =0.347+ 0.010°K for T less than 0.1°K, exactly 














i} T T T 7 
q Runt S | 
4 . Run2 9 
a Oo o Run3 o » 
. Run4 V7 
° o o 
lo © = al 
o.)6|U° ° ° Vv o 
° 9 o V Vv 
l a | l | 
10 18 26 34 42 50 








VoLuME 7, NUMBER 6 


PHYSICAL REVIEW LETTERS 


SEPTEMBER 15, 1961 





the result given by AHW. Within each run we 
found x to be constant to within about 2%, so it 
is likely that x is more nearly constant than is 
indicated by all the data collected in Fig. 2, for 
which there is about 5% scatter. The slight drop 
in x below 0.05°K observed by AHW and attributed 
to experimental effects was not observed in these 
measurements. 

The susceptibility can be written [Eq. (5.8), ref- 
erence 5] 


B? Ss (ar /€),, p2? N 


x*7Tl +TE(@r/ee) 4 VaT* (1) 


an equation which defines T*. In this equation N/V 
is the number of He*® atoms per unit volume, B=, 
is the product of the gyromagnetic ratio and 
Planck’s constant divided by 21, k is Boltzmann’s 
constant, (87/8€), is the density of states at the 
Fermi surface, and {is the average of the effective 
spin interaction function over the Fermi surface. 
Using the most recent empirical values for the mo- 
lar volume”?! and effective mass,’ the value of 
T* given above leads to [1 + 4£(87 /8€),, ]= 0.293. 
Hone’ writes the diffusion coefficient in the form 


D= $V) Tpit + 10(07/d¢) L (2) 


where V,=P,/m* is the quasi-particle velocity at 
the Fermi surface, Tp is a relaxation time for 
diffusion, and the quantity in square brackets is 
given above. The resulting empirical value of Tp 
is 4.6x107°T~ sec °K?. If one uses in addition 
to the above information the most recent value 
for the speed of sound" to estimate’ the value of 
Tp from the Fermi liquid theory, one obtains tp 
=1.0x107*°T™ sec °K*. A relaxation time can 
also be computed from the measured values of 
thermal conductivity® in the Fermi liquid region 
giving Tx =6.4x107"* T™ sec °K*. If there were 
no dependence of the quasi-particle scattering 
probability on the azimuthal angle of scattering 
and in addition no spin effects, then by theory’’’ 
one has tp/tx = 4; experimentally we observe 

Tp /t K = 9.70. 
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Using the above values of the relaxation times, 
one finds in both cases that the ratio of the uncer- 
tainty in quasi-particle energy to kT is less than 
1 at 0.04°K: #/kTr =0.5 and 0.7 using Tx and Tp, 
respectively. Moreover, since the specific heat, 
thermal conductivity, and self-diffusion coefficient 
all have the predicted temperature dependences be- 
low 0.04°K, we conclude that the concept of He’ as 
a Fermi liquid at temperatures less than this and 
at low pressures is valid. 

We wish to acknowledge helpful discussions with 
Professor John Bardeen and D. Hone. We would 
like to thank William Abel for his expert assist- 
ance in constructing apparatus and for his help 
in making the measurements. 
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OBSERVATIONS OF IONIC SOUND WAVES IN PLASMAS 


I. Alexeff and R. V. Neidigh 
Oak Ridge National Laboratory, Oak Ridge, Tennessee 
(Received July 17, 1961; revised manuscript received August 31, 1961) 


The object of this Letter is to report the ob- 
servation of ionic sound waves in a typical dis- 
charge-tube kind of plasma, and also in a plasma 
restricted by the presence of a magnetic field. 
Earlier, partial reports of these experiments 
have been given local distribution.** The earliest 
work in this field is thought to be the prediction 
of ionic sound waves in plasmas by Tonks and 
Langmuir in 1929. The fundamental frequency 
and first and second overtones were found in a 
mercury vapor discharge confined in a spherical 
tube by Revans in 1933.5 He also showed evidence 
for radial standing waves in a cylindrical tube. 
Data from the recent discharge-tube experiments 
of Crawford® have been analyzed by this labora- 
tory' and by Moore’ and shown apparently to ex- 
hibit ionic sound waves. Additional work by Craw- 
ford® confirms this analysis. 

The theory of propagation of ionic sound waves 
is much like that of ordinary sound waves. The 
velocity v of an ionic sound wave is related to the 
electron temperature FT, and the ion mass m; by 
the expression 


v= (kT, /m)”. 


Here k is Boltzmann’s constant and y is the adia- 
batic compression coefficient which is § for or- 
dinary sound waves and may be as high as 3 in 
rare plasmas.® 

If the plasma is restricted to a long cylindrical 
column by a magnetic field and the sound wave is 
reflected similarly at each end, standing ionic 
sound waves may result as ordinary sound waves 
do in an open organ pipe. Then the fundamental 
frequency should be v/2L, where L is the length 
of the column, and the harmonics should be inte- 
gral multiples of the fundamental. 

If the plasma is confined by a spherical dis- 
charge tube without a magnetic field, the fre- 
quency of the fundamental is v/CD, where D is 
the sphere diameter, and C is a constant that 
can be 1.51 or 1.0 depending on the plasma 
boundary conditions. Higher frequency modes 
of oscillation should occur, but their frequency 
also depends on the boundary conditions. Fora 
pressure antinode at the wall (C =1.51), the fre- 
quency ratios for the first six overtones to the 
fundamental are 1.61, 2.16, 2.17, 2.71, 2.85, 
3.25. 


A schematic diagram of the apparatus which 
produced the cylindrical plasma column is shown 
in Fig. 1 (inset). The plasma column was pro- 
duced by a }-inch diameter electron stream (up 
to one ampere) from a negatively biased (up to 
-120 volts) hot filament. The magnetic field 
could be varied from 1000 to 7000 oersteds. The 
frequency as observed on a wire probing the col- 
umn was found to be grossly independent of the 
magnetic field and the filament bias. It did not 
depend on background gas pressure which was 
uniform over the length of the column and in the 
range 10~ to 10-° mm Hg. The filament temper- 
ature had the most effect on the frequency which 
could be “pulled” slightly but which would change 
discontinuously to a higher mode if “pulled” too 
far. The amplitude of the probe signal varied 
considerably (more than 20 volts) with change 
in the above parameters. Sinusoidal frequencies 
were observed up to the fifth harmonic. At times 
the wave shape gave evidence of higher harmon- 
ics present with the fundamental. 

The observed variation of frequency with ion 
mass for two lengths of plasma column is shown 
in Fig. 1. The solid lines are the theoretical 
curves for ay of 3, and an electron temperature 
of 17 ev, the average of the first excitation po- 
tential of the inert gases. Electron temperature 
measurements made with a Langmuir probe were 
within a factor of two of this value. Obtaining re- 
liable probe characteristic curves proved to be 
very difficult under the conditions of operation. 
In general, gases whose first excitation potential 
was greater than 17 ev supported standing waves 
with frequencies higher than predicted, and lower 
excitation potentials resulted in correspondingly 
lower frequencies. 

For a convincing proof of the existence of these 
waves it was felt that they should be exhibited in 
a different system, one not using a magnetic field. 
After preliminary experiments, a spherical glass 
discharge tube was chosen instead of the usual 
cylindrical type. It was thought that the lowest 
mode of oscillation in a sphere might be less 
damped, as this mode involves no transport of 
matter parallel to the wall. The diagram of the 
sphere is shown schematically in Fig. 2 (inset). 
The sphere is glass with a hot cathode of tungsten 
wire and a metal disk as an anode. Both are flush 
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ATOMIC MASS 


FIG, 1. Frequency vs atomic 
mass for two lengths of plasma 
column, The inset is a sche- 
matic diagram of the 23-cm arc, 
The 2,2-cm arc was excited in 
the same manner. 


FIG, 2. Frequency vs atomic 
mass for two diameters of 
spheres, The inset is a sche- 
matic diagram of the sphere 
geometry. Electrode sizes 
were proportional to sphere 
diameter in each size sphere. 
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Table I. Ratios of frequencies observed to the fundamental for the different gases used. The calculated ratios 


are given in the left-hand column. The lowest frequency 


observed for helium, neon, and argon was thought to be 


the first overtone since this provided a better fit of the data. 











C=1.51 Hydrogen Helium Neon Argon Krypton Xenon 
1 1 
1.61 1.46 1 #1 1.61 1.61 1.75 1.57 
2.16 1.93 2.00 2.16 
2.17 2.30 2.22 2.20 
2.71 2.72 2.59 2.63 
2.85 
3.25 3.43 








with the wall and on opposite ends of a diameter. 
A Langmuir probe in the center was used to meas- 
ure the electron temperature. Oscillations were 
detected by means of a capacitively coupled elec- 
trode outside of the sphere. Pressures of about 
10-? mm Hg and currents of 100 ma resulted in 
oscillation amplitudes of about one volt in the de- 
tection circuit. Varying the discharge param- 
eters caused the oscillation frequency to change 
discontinuously to the higher frequency modes. 
The lowest frequency mode that was thought to 
correspond to the fundamental of the sphere was 
very difficult and sometimes impossible to ex- 
cite. Probably some asymmetry in the sphere 
preferentially excited higher frequencies, and 
the high gas pressure needed to maintain the 
discharge in the sphere preferentially damped 

the lower frequencies. 

In analogy with the plasma column work the low- 
est frequency observed in plasmas of different 
atomic mass and for spheres of 2.5 and 25 cm 
in diameter are plotted in Fig. 2. As the meas- 
ured electron temperature now varied widely, all 
frequencies were normalized to an electron tem- 
perature of 1 ev. The solid lines represent the 
lines of slope (-4$) which pass through the lowest 
observed frequency for both sphere diameters. 
Assuming that the lines pass through the funda- 
mental frequencies, y = § and C=1.51. The value 
y= is reasonable for the spheres in view of the 
high gas pressure in the discharge. In general 
the frequency is seen to be proportional to (D)™ 
and (m;)““*. Apparently the fundamental was not 
excited in some of the gases. 

As many as four different frequencies have been 
observed for one gas. Classification is difficult 
because the fundamental might not have been ex- 
cited and the frequency on an overtone may be 


altered by possible changes in the boundary con- 
ditions. An attempt at classification appears in 
Table I. 

Another complication in working with ionic 
sound waves in discharge tubes appears if the 
pressure is too high. At about 1 mm Hg strong 
oscillations appear which are lower in frequency 
than the fundamental of the ionic sound waves. 
These frequencies apparently correspond to or- 
dinary sound waves in the residual gas. 

The authors gratefully acknowledge the sugges - 
tions, impetus, and continued support given the 
program by A. H. Snell and E. D. Shipley and the 
valuable discussions with many others of the 
Thermonuclear Division staff. 
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SUPERCONDUCTING-NORMAL TRANSITION TIME 


A. H. Nethercot, Jr. 


Thomas J. Watson Research Center, International Business Machines Corporation, Yorktown Heights, New York 
(Received August 9, 1961) 


The nonequilibrium properties of superconduc- 
tors are of particular interest since these prop- 
erties cannot be calculated directly from the 
standard Bardeen-Cooper -Schrieffer (BCS) the- 
ory. One such problem is the calculation of the 
speed at which transitions can occur between 
the normal and superconducting states when the 
magnetic field is varied. Since the theory does 
not clearly indicate what fundamental effects 
might limit this speed, a number of possible ex- 
planations have been suggested between which it 
is difficult to choose (see Appendix). It is the 
purpose of this Letter to propose a new explana- 
tion which proceeds in a relatively direct and 
straightforward fashion from the BCS theory 
and which, since it predicts a rather long time, 
overrides many of the previous suggestions. 

In brief, the present explanation is that the 
switching speed is limited by a spatial rather 
than by a purely temporal effect: At fast switch- 
ing speeds and therefore at high frequencies the 
skin depth becomes so small that this thin layer 
is prevented from changing its state by the pres- 
ence of the large amounts of unswitched adjacent 
material. As will be described below, it appears 
from present theories’? that skin depths corre- 
sponding to switching speeds of 10~*° sec are so 
thin that their properties are severely altered 
by the presence of the adjacent unswitched ma- 
terial. Hence these skin depths cannot be said 
to change their state sufficiently to define switch- 
ing in times on the order of 107*° sec. 

The predictions of this theory are in accord 
with the presently known experimental facts. It 
is known from experiments with pulses that super- 
conductors do switch in at least several millimi- 
croseconds. However, experiments*” on the non- 
linear properties of superconducting tin at 10000 
Mc/sec indicate that switching is only very in- 
complete or partial in times on the order of 2 
x10 sec. Thus the switching time is bracketed 
between 10~° and 107" sec. In these 10000-Mc/ 
sec experiments, the amount of second harmonic 
power generated when a high-amplitude rf mag- 
netic field is incident against a bulk or thin film 
superconductor is compared with the harmonic 
power calculated under the assumption that there 
is an instantaneous change in the microwave sur- 
face impedance from its normal state value to its 
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superconducting value when the field changes 
from above critical to below critical during a 
cycle. Distortion of the reflected wave is ob- 
served only when the rf field crosses the critical 
field curve during a cycle. The amount of this 
distortion is less than one percent of that calcu- 
lated. This indicates that there is only a partial 
change in state. In addition, below 3.2°K harmon- 
ics are observed only if the bias field, H,, is 
greater than H,, indicating that at low tempera- 
tures switching from superconducting to normal 
is much more difficult than switching from normal 
to superconducting. Also, the amount of harmon- 
ics observed from thin films decreases very rap- 
idly as the films become thinner than 2500 A. 
Thus, such films do not switch as readily as bulk 
materials, at least when the switching field is ap- 
plied to only one side of the film. 

The theoretical situation can be discussed from 
either the point of view of Parmenter’ or that of 
Cooper.” The former shows that for supercon- 
ducting-normal-superconducting (S-N-S) or N-S-N 
sandwiches there exists a critical distance, 2X;, 
for the central portion. For tin 2X-=1500 A for 
both the S-N-S and N-S-N sandwiches. For thick- 
nesses, d, smaller than 2X; the central region 
shows the general properties only of the adjacent 
material, while for thicknesses, d, greater than 
2X; the central region displays its own properties 
(but if it is superconducting it has a reduced ener- 
gy gap even if d>2X;,). We will discuss first the 
case in which a bias field, H,, is applied to a 
semi-infinite bulk superconductor, where H9>H-- 
An additional microwave field, H,¢, of large am- 
plitude is also applied to the surface. The effect 
of this field is to tend to restore the skin depth to 
the superconducting state during one-half of the 
cycle. Therefore we wish to consider a thin layer 
one skin depth thick of superconducting metal in 
contact with a bulk normal metal (if 0 denotes 
empty space, an N-S-0 sandwich). Although the 
N-S-0 sandwich is not treated by Parmenter, it 
would be expected that the central region could 
develop a gap only if the thickness is greater than 
a critical length, 2X,’. This distance, 2X¢’, 
would be expected to be somewhat less than 2X, 
since the presence of N material on only one side 
rather than on both would be less effective in pre- 
venting formation of the gap. A rough estimate 
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leads to a value for 2X,’ of 0.8(2X,,) or 1200 A 
for the N-S-0 sandwich. If the thickness is only 
slightly greater than 2X’, the average gap in the 
central region is =1/30 bulk gap.' Any gap, no 
matter how small, suffices to define the super- 
conducting state for the dc currents. However, 

at high frequencies and T>0, the gap, AE, must 
be greater than a certain amount for some prop- 
erty such as the microwave surface resistance 

to decrease by a factor of 2 from its normal 

state value. From inspection of data on the mi- 
crowave surface resistance versus temperature, 
this decrease of a factor 2 is seen to occur when 
kT=AE. If we consider the temperature range 
from 3.7 to 2.5°K, it follows that the gap must be 
less than the bulk gap only by a factor five (rather 
than by a factor thirty) for this given decrease in 
surface resistance to occur. This large a gap 
will occur only at distances rather greater than 
2X’, estimated to be 3 to 4X,.’.. Thus the normal 
skin depth (5,7) must be greater than 2000 A if 
this depth is to be able to change its state from 
NtoS. The skin depth becomes approximately 
2000 A at v=3x10° cps, and hence it is predicted 
that switching cannot occur when the field is 
changed from greater than critical to less than 
critical in times shorter than 7 = 1/2my = 0.5 x 107° 
sec. 

It is predicted that switching from the S to N 
state is even slower than the above. When Hy<H ° 
the penetration of the field is governed by the su- 
perconducting skin depth, 55. Hence, at the lower 
temperatures where 5¢=A, switching will be ex- 
tremely difficult since A (the penetration depth) 
is only 500 A. This is in agreement with experi- 
ment* since no nonlinear effects are seen in tin 
at 10000 Mc/sec if Hg<H, and T<3.2°. It might 
be objected that this approach leads at low tem- 
peratures to a frequency-independent limitation 
on switching. However, at sufficiently low fre- 
quencies new effects can occur such as motion 
of S-N boundaries and nucleation and growth of 
new sites. These effects will then allow switch- 
ing to take place. 

In the above, the S-N-0 case has been differen- 
tiated from the N-S-0 case with the important 
parameters being 5s and 65y, respectively. If 
switching is complete, the spatial variation in 
both N and S material should be included in some 
complex manner. However, as a first approxima- 
tion, a high enough frequency can be chosen that 
Switching is only partial and hence only one skin 
depth is important. 

The above treatment can also be carried through 


on the basis of Cooper’s theory.” This theory does 
not give a critical distance below which the gap 
completely disappears in a superconducting layer, 
but instead the gap decreases in size indefinitely 
as the layer becomes thinner. Adopting again the 
criterion that the gap cannot be less than one-fifth 
of its bulk value (to establish a change of a factor 
two in the microwave surface resistance) and dis- 
cussing the case Hy>H, (the N-S-0 sandwich), 

we find from Eqs. (1) and (3) of reference 2 that 


exp[-N(0)Vt,,/t,] =4, 


where ts is the skin depth and ty is the effective 
thickness of the normal material. Since N(0)V 

= 0.30 for tin, tn /ts =4. If a limit to (ty +ts) is 
set by the bulk correlation distance (&,=2500 A), 
we find that ts > 1700 A and that the switching 
time, 7, is longer than 0.3x10~° sec. This value 
is in reasonable agreement with the 0.5x107*° 
sec from Parmenter’s theory. However, it ap- 
pears that this value may be too short since, from 
the uncertainty principle, the correlation distance 
for reasonably pure bulk materials must be in- 
versely proportional to the energy gap. If the 
modified gap is } its bulk value, the correlation 
distance is 12000 A and this leads to a switching 
speed of 3x10~™ sec, a very slow speed. Cooper’s 
theory needs more development to clarify which 
correlation distance is involved. 

The theory of Cooper is particularly suited for 
the case of thin films. If the switching field is 
applied to only one side of the film and if the film 
is thinner than the skin depth (either 55 or 5y), 
the microwave field varies approximately linearly 
across the film and falls to a very small value on 
the other side. Therefore, in general, the micro- 
wave field falls off more rapidly than it does in 
bulk material and at a given frequency the inter- 
action potential N(0)V is more diluted and the gap 
is less. Thus the degree of switching would be 
less in the thin film. These results for thin films 
are in accord with experiment* since nonlinear 
effects at v=10000 Mc/sec rapidly decrease in 
amplitude for thicknesses less than =2500 A. 

There is no experimental evidence on the case 
when the microwave field is applied to both sides 
of the thin film. If Hg<H,’, where H,’ is the thin- 
film critical field, the situation is very similar 
to the above except that the field varies even more 
rapidly in space. The switching should be very 
slow. The other case, Hg>H,’, is more difficult 
to treat since the fields do not vary greatly in 
space. However, the Gor’kov theory indicates 
that the transition becomes second order in the 
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magnetic field for the thinner films® and this 
causes a limitation on the switching. For first 
order transitions an infinitesimal change of H, 
from (H,.’+€) to (H,’- €) causes a well-developed 
gap to appear and a large change in the micro- 
wave surface resistance to occur. In a second 
order transition a much larger change in field is 
necessary to produce the requisite gap. It is 
clearly much more difficult to produce a large 
change in field in a short time (=107*° sec) than 
it is an infinitesimal change. The switching field 
should be at least several percent of the critical 
field to produce a gap } of the bulk gap.* Since 
the critical field for thin films is very much 
greater than that for bulk materials, the requisite 
field may thus be quite appreciable. A clearer 
understanding of the relation between the Gor’kov 
and the spatial dependence theories would lead to 
a better understanding of these temporal effects 
as well as of the high-field superconductivity ef- 
fects associated with filamentary behavior. 

We have had informative and interesting conver - 
sations with a number of people. Discussions 
with R. J. von Gutfeld and R. H. Parmenter have 
been particularly helpful. 

Appendix. Since no discussion of the various 
possible explanations for a limitation of the switch- 
ing time has been given previously, a brief listing 
of some of these will be given here. The usual lim- 
itation, that of the growth of normal domains lim- 
ited by eddy currents, has no relevance when skin 
depths small compared to the correlation distance 
are involved. 

(I) The uncertainty relation AEAt=f leads to 
At=107™ sec if AE is the energy gap. This time 
is too short. 

(II) The correlation distance divided by the pho- 
non velocity (£9/v,) leads to a time of about 10~*° 
sec and could be interpreted as a signaling time 
between correlated electrons by means of the in- 
teraction responsible for superconductivity. How- 
ever, it is difficult to guess from the BCS theory 
whether such an effect should in principle be 
present and in addition the electrons are physi- 
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cally moving much faster than the phonons. 

(III) A relaxation time such as l/vg (where / is 
the electron mean free path and uf the Fermi 
velocity) should come into the theory since an 
electron trajectory cannot be greatly altered un- 
til it is terminated. However, in the treatment 
of the anomalous skin effect by the ineffectiveness 
concept it is not always clear whether //vf or 6/v 
should be used.® The former would lead to switch- 
ing times of from 10~° sec to much shorter values 
and the latter to about 107"* sec. The experimen- 
tal fact that similar switching behavior is ob- 
served in situations where / should be drastically 
different indicates that L/vf may not come in di- 
rectly. 

(IV) It has also been speculated’ that the velocity 
of second sound (vqq) in the electron ensemble 
could possibly limit the rate of conversion of su- 
perconducting electrons to normal. This would 
lead to a time 7 =A/vyz=10™ sec if* vyy=[H,70 -f*/ 
27p}=10* cm/sec. However, it is not clear that 
relaxation effects would allow this type of second 
sound to be propagated (especially in impure met- 
als or in evaporated films) or even whether sec- 
ond sound of a different variety® and with a veloc- 
ity of ~10° cm/sec would or would not be pre- 
ferred. 





'R. H. Parmenter, Phys. Rev. 118, 1173 (1960). 
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3A, H. Nethercot, Jr., Proceedings of the Seventh 
International Conference on Low-Temperature Physics 
(University of Toronto Press, Toronto, 1960), p. 231. 
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5D. H. Douglass, Jr., Phys. Rev. Letters 7, 14 
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8See J. Bardeen and J. R. Schrieffer, in Progress in 
Low-Temperature Physics (North-Holland Publishing 
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*D, J. Thouless and D. R. Tilley, Proc. Phys. Soc. 
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TWO-PHOTON EXCITATION IN CaF,:Eu** 


W. Kaiser and C. G. B. Garrett 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received August 28, 1961) 


With the development of optical masers,!~° it 
is now possible to study two-photon processes 
which necessitate intense sources of monochro- 
matic radiation. We have investigated the genera- 
tion of blue fluorescent light around \,= 4250 A 
by illuminating CaF,:Eu** crystals with red light, 
\, = 6943 A, of a ruby optical maser.” Our experi- 
ments differ from the recent investigations by 
Franken et al.,* who observed the generation of 
optical harmonics in quartz. It is essential in 
their experiments that the crystal lacks a center 
of inversion and that it is transparent at Vy and 
2v,. In our investigations, the highly symmetric 
cubic CaF, structure is used and after an excita- 
tion to a real absorbing state of Eu** at 2v,, the 
fluorescent decay to a lower state is observed. 
CaF, crystals, with 0.1% Eu®* ions substituted 
for Ca**, exhibit strong absorption between 30000 
and 25000 cm™ resulting from electronic 4f + 5d 
transitions of the Eu** ion.’ Excitation of 
CaF,:Eu** with light absorbed in this wavelength 
range gives rise to a brilliant blue fluorescence 
(bandwidth ~ 300 A) which originates at levels 
around 4200 A and terminates at the ground state. 
There are no energy levels below 22000 cm™, 
making CaF,:Eu®* crystals transparent at v= v,. 

We have performed the following experiments. 
A CaF,:Eu* crystal of 1-mm thickness was 
mounted in front of the entrance slit of a quartz 
spectrometer, and the light beam of a ruby opti- 
cal maser (~0.1-joule output) was focused onto 
the sample. Two red filters with transmission 
values T<10"* for , <6100 A were mounted be- 
tween the maser and the CaF,:Eu** crystal in 
order to eliminate the possibility that blue or ul- 
traviolet stray radiation from the flash lamp 
reaches the sample. The fluorescence of the 
crystal was photographed at the exit slit of the 
spectrometer. Between the sample and the prism, 
a blue filter with T=10~ at 7000 A was interposed 
to avoid excessive blackening of the photographic 
plate. In Fig. 1 a picture of the output spectrum 
is presented which was obtained by one flash 
(~500 sec) of red light. The overexposed bright 
Spot at 6943 A results from the incident maser 
light. Of particular interest is the observation 
of light around 4250 A, which is characteristic 
for the blue fluorescence of CaF,:Eu**. The width 
of the blue streak is approximately 0.3 mm cor- 


responding to the diameter of the incident light 
beam. When pure CaF, was illuminated by the 
optical maser in the same way, no light with 

\ < A, Was observed on the photographic plate. 
This observation is expected from the high sym- 
metry of the CaF, lattice. 

In a second experiment, the light beam leaving 
the quartz spectrometer was intersected by two 
mirrors in such a way that the red and blue parts 
of the spectrum were directed separately onto 
two photomultipliers. The signals of these photo- 
multipliers were displayed simultaneously on a 
dual-beam oscilloscope. Load resistors of 10° 
ohms in the photomultiplier circuit were employed 
in order to smooth out the spikes (resulting from 
the relaxation oscillations) of the ruby optical 
maser and to allow a direct quantitative compari- 
son between the two signals. In Fig. 2 the signal 
I; obtained in the blue part of the spectrum is 
plotted against the signal in the red J,, which is 
a direct measure of intensity incident on the crys- 
tal. The empirical line through our experimental 
points represents the quadratic relation J b «l?, 
which is a strong indication that we are dealing 
with a two-photon process. From the calibration 
of our photomultiplier, the number of blue photons 
at the exit slit of the spectrometer could be esti- 
mated to be 5x10’ photons per flash. This number 
agrees well with the value of ~10® photons per 
flash estimated from the blackening of the photo- 
graphic plate shown in Fig. 1. 

We believe that in our experiments a real state 
of the Eu** ion at v=2v, is excited by a two-photon 
process. The observed fluorescence is only an in- 
dication that the ion was indeed in the excited 
state. The probability of a process by which an 
atom is excited by two photons hy, andhv,, the 
sum of which corresponds to an excited state of 
the atom, was first treated by Goeppert-Mayer in 


5000 6000 8000 


3000 4000 





FIG. 1. Positive of photographic plate, indicating 
the blue emission of a CaF,:Eu* crystal under strong 
illumination with A,=6943 A. 
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FIG. 2. Blue fluorescent in- 
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1931.° More recently Kleinman’ has reinvesti- 
gated the problem and has derived an equation for 
the cross section o for such a two-photon excita- 
tion process which allows a direct quantitative es- 
timate of the expected number of blue photons: 


o= (e?/me*Pa?/n*a v)F, (1) 


where e?/mc*=2.8x107"* cm, n is the dielectric 
constant of the material, Av is the width of the 
absorption band around 2v,, and F is the incident 
photon flux. The number of excitations or, as- 
suming a fluorescent quantum efficiency of one, 
the number of emitted fluorescent photons per 
unit volume of the specimen is 


P a" ONF = e?/me*?Q,2/n*a v)NF?. (2) 


Here N denotes the number of fluorescent ions in 
230 


the crystal. For an estimate of P», the following 
experimental data were used: ,=7x10~* cm, 
n(CaF,)=1.4, Av=1.5x10"* sec™ (5x10° cm™), 
and N= 2.4x10"* Eu ions per cm*. The maser out- 
put of 4x10*” photons was concentrated into an 
area of approximately 107° cm?; the incident pho- 
ton flux F is, therefore, F = 8x10** photon cm™ 
sec™, From Eq. (2) we obtain P,=10** photons 
cm™ sec™. Since the active volume was 10~™ 
cm’, one calculates 5x10*° blue photons to be 
emitted per flash (500 usec). The efficiency of 
this two-photon process is therefore ~10~. This 
value could be greatly increased if a system with 
a narrower excited state (smaller Av) were em- 
ployed. For a comparison with the observed num- 
ber of blue photons, one has to consider the f/4 
lens system and a number of reflection losses in 
the spectrometer. The number of photons leaving 
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the instrument is calculated to be 10°. This value 
agrees quite well with the experimentally observed 
photon flux, giving strong support to the interpre- 
tation of our observations. The agreement is bet- 
ter than expected considering the simplifications 
made in the derivation of Eq. (2) and the estimates 
made during its application. 

We would like to thank D. Kleinman for commu- 
nication of his paper prior to publication, and 
H. Guggenheim for growing the CaF,:Eu** crystal. 
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MAGNETIC BREAKDOWN IN CRYSTALS* 


Morrel H. Cohen and L. M. Falicov 


Institute for the Study of Metals, University of Chicago, Chicago, Illinois 
(Received June 7, 1961) 


To establish the various regimes of behavior of 
electrons in a magnetic field, the cyclotron fre- 
quency w, is usually compared with tr“, A“kT, 
ork“ Ep, where 7 is a typical relaxation time 
and Ey the Fermi energy. Comparison of w, 
with AE g» where Eg is an energy gap of the 
pertinent band structure, has not been discussed 
except to note that the usual derivations of an ef- 
fective Hamiltonian’ for the motion of electrons 
in a magnetic field is no longer valid when iw, 
~Eg. No attention has been directed to the ques- 
tion of what actually happens when fiw, exceeds 
E gs that is the subject of the present Letter. 

The splitting of the energy levels due to mag- 
netic fields is of the order of 10~ ev per gauss, 
whereas gaps in ordinary band structures are 
rarely as small as 0.1 ev. Our question would 
be academic were it not for the existence of small 
gaps of the order 107° ev or less in some metals 
because of spin-orbit splitting? or of points of ac- 
cidental degeneracy at or near the Fermi level.° 
These gapd ‘are small compared with the Fermi 
energy and the remaining typical gaps in the band 
structure. In such a case the answer to our ques- 
tion is that as H increases, one passes from elec- 
tron orbits determined semiclassically from the 
entire band structure to those determined semi- 
Classically by ignoring the small gaps. We call 
this effect magnetic breakdown. Of course, some 
additional effects take place, i.e., a spreading of 
the magnetic levels into narrow bands and addi- 
tional scattering processes, but these are not sig- 
nificant except in the transition range of magnetic 
field Rw, ~Eg- 





Let us consider a simple example of a free- 
electron metal perturbed by a one-dimensional 
periodic potential. The one-electron Hamiltonian 
in the presence of a magnetic field in the z direc- 
tion is 


1 lel+ 
Kem 5 {B+ A) + Vo cosne, 
where 
A=(0, Hx, 0), Re >K. 


Since the electron distribution overlaps into the 
second zone, the Fermi surface consists of an 
undulating cylinder in the first zone with axis 

in the x direction and a small pocket of electrons 
in the second zone [Fig. 1(b)]. When H is small, 
i.e., iw-«Vo, the undulating cylinder gives rise 
to a set of open orbits; as consequences de Haas— 
van Alphen oscillations arise from the pocket but 
not from the cylinder and the R,, component of 
the resistivity tensor increases quadratically with 
H. When hw,.>Vo,* magnetic breakdown essen- 
tially restores the free-electron surface [Fig. 1(a)] 
and only the equatorial section of the “sphere” 
gives rise to a de Haas—van Alphen period. More- 
over, because the open orbits have now disap- 
peared, R,, decreases as H™~ to a saturation 
value which is very approximately the free-elec- 
tron value. The spreading of the magnetic levels 
into bands, which have a maximum width of about 
2V, for orbits which just intersect the zone faces, 
gives no appreciable effect, and the additional 
scattering processes give only second order cor- 
rections to the free-electron magnetoresistance® 
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A simple model to illustrate magnetic break- 


FIG, 1. 
down, (a) shows the orbits in k space for the free-elec- 
tron case which are the same as the orbits when mag- 
netic breakdown occurs; (b) shows the orbits for the low- 
field case, 


that vanish as H~™. 

We have seen that in this simple illustration 
magnetic breakdown changes open orbits into 
closed orbits and so the H? term in the magneto- 
resistance disappears. More generally, the ef- 
fect can cause a change of any orbit, i.e., closed, 
extended, or open, into any other of the same or 
different category. Consequently, as the field 
varies it must be possible to observe the disap- 
pearance of some periods and the appearance of 
new ones in experiments such as the de Haas— 
van Alphen effect, cyclotron resonance, and ul- 
trasonic attenuation. 

Possible experimental evidence for magnetic 
breakdown is the giant orbit observed by Priest- 
ley® in the de Haas—van Alphen effect in magne- 
sium for fields very near the (000, 1] direction. 
Such an orbit, with an area larger than the sec- 
tion of the Brillouin zone, cannot occur in any 
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ordinary band structure. However, the existence 
of an accidental degeneracy near the Fermi ener- 
gy along the IK line (without spin-orbit coupling) 
and the detailed analysis of the structure of the 
levels*® shows that the waist of the Fermi surface 
in the second zone contacts each “cigar” in the 
third zone just above and below the basal plane. 
Spin-orbit coupling should lift the degeneracy, 

but a sizeable region of k space must remain in 
which magnetic breakdown occurs for orbits which 
pass near one set of six points of contact. This 
would cause transitions between the Fermi surface 
in the second zone and the “cigars” in the third 
zone, giving rise to an orbit essentially the same 
as the equatorial section of the free-electron 
sphere (Fig. 2); its area gives a value of the 
period that exactly agrees with the experimental 
value 1.57x10~ (gauss)~*. The narrow angular 
region (about 4° in radius around the c axis) in 
which the orbit has been found further supports 


























FIG, 2. The giant orbit in magnesium, The figure 
shows a section of the Fermi surface with a plane paral- 
lel to the basal plane and passing through the points of 


contact removed by spin-orbit splitting. --— orbit 
ardund the Fermi surface in the second zone at low 
fields; +++ orbit around “cigars” in the third zone at 
low fields; —~ giant orbit after magnetic breakdown oc- 
curs, showing transitions from second zone to third 
zone to second zone, etc. 
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this interpretation, because when the field is 
tipped only slightly off the c axis, the orbit can 
no longer pass near all six points of contact. The 
fields used in the experiments (57 to 191 kgauss) 
should have been enough to overcome the small 
gaps, but we have not estimated the splitting of 
the accidental degeneracy. 

Another case of magnetic breakdown which may 
occur at relatively low fields is the magnetore- 
sistance of beryllium for fields perpendicular 
to the c axis. If the H, level’ of the band struc- 
ture lies below the Fermi energy, spin-orbit 
coupling? causes the Fermi surface in the third 
zone to look like two similar cylinders with axes 
parallel to the c axis; the magnetoresistance 
must then increase quadratically with H. How- 
ever, for fields of about 20 kgauss, magnetic 
breakdown connects these cylinders with the 
small pockets in the fourth zone, overcoming 
the small energy gap produced by spin-orbit 
splitting. The orbits close, and the magneto- 
resistance must saturate. A similar effect should 
occur in magnesium at fields very roughly about 
200 kgauss from the change in topology of the 
Fermi surface in the second zone.? 

In Mg,® Zn,® and Cd® there are two intersecting 
electron pieces of Fermi surface around L, which 
are separated by spin-orbit coupling, except at 
two points of contact on the AL line. For the mag- 
netic field in a plane perpendicular to AL, mag- 
netic breakdown is possible for orbits near the 
points of contact at any field strength. For arbi- 
trary field orientation hw, must exceed the spin- 
orbit splitting. In the hexagonal close-packed 
phases of Li and Na, the piece of Fermi surface 
around A which is split off from the main body 
of the Fermi sphere by spin-orbit coupling must 
show similar behavior because of the vanishing 
of the splitting along the three lines LAL. 

The points of contact on the Fermi surface of 
graphite® are also removed by spin-orbit split- 
ting.*° Magnetic breakdown phenomena should 
in principle be observable when the field is nearly 
normal to the c axis. 

Finally, it should be remarked that band gaps 
introduced by spin-orbit coupling can be magnetic 
field dependent via the interaction of the magnetic 


field with the unquenched orbital angular momen- 
tum and the electron spin. For the simple exam- 
ple of the degeneracy of p-like states in a reflec- 
tion plane, the condition for magnetic breakdown 
is unchanged when the magnetic field is in the 
plane. On the other hand, magnetic breakdown 
can occur at somewhat different field strengths 
for the two different spin orientations when the 
magnetic field is perpendicular to the plane. 

In conclusion we would like to acknowledge 
several stimulating discussions on this subject 
with W. A. Harrison. 





*Work supported in part by the Office of Naval Re- 
search. 
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MULTICOMPONENT MAGNETOPLASMA RESONANCE IN p-TYPE GERMANIUM 


R. E. Michel and B. Rosenblum 
RCA Laboratories, Princeton, New Jersey 
(Received July 31, 1961) 


In a field of sufficiently high frequency w the 
electrical conductivity of a system of carriers 
exhibits a maximum (cyclotron resonance) when 
we =eH,/m*c=w, where H, is the applied mag- 
netic field and m”* is the effective mass. At car- 
rier concentrations high enough that the plasma 
frequency of the carriers is comparable to the 
experimental frequency, the usual cyclotron ab- 
sorption maxima are shifted by plasma effects. 
We have observed such phenomena in p-type 
germanium where the plasmas of the light and 
heavy holes form a coupled system. The study 
of such a coupled system is particularly inter- 
esting in that it gives insight into certain aspects 
of carrier motion in metals and semimetals. 

Magnetoplasma resonance was first observed 
in solids in the case of indium antimonide with 
an approximately fixed concentration of a single 
type of carrier.’ It has also been studied in n- 
type germanium’ as a function of carrier concen- 
tration. 

In the present experiments a small, approxi- 
mately ellipsoidal, sample of 7 ohm-cm indium- 
doped germanium was placed in a microwave cavi- 
ty and its conductivity at 24 kMc/sec determined 
as a function of magnetic field in the [110] direc- 
tion. The number of holes in the valence band was 
controlled by adjusting the temperature in the 10 
to 40°K region. The carrier concentration was 
determined by Hall effect measurements. The el- 
lipsoid for the plasma experiments was cut from 
between the potential arms of the Hall sample to 
minimize error due to doping inhomogeneities. 

In an ellipsoidal sample which is small com- 
pared to the wavelength of the applied field and 
also small compared to the penetration depth, the 
microwave electric field will be uniform over the 
entire sample. For this case a simple equation 
of motion for the carriers has been derived,’ Kit- 
tel® has extended this treatment to the case of a 
plasma consisting of two mobile particles and the 
fixed ions. The two coupled equations of motion 
can be written: 


= 1: ; 2< mM, 2+ _¢s .% + 
1s)" ae ta. tet 
a . 34 2 m 2 ez > = 
—_— +—ty =e 
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where fr, and fr, are the position vectors for the 
average particle of the first and second kind (in 
the present case the light and heavy holes, re- 
spectively); m, 2 are their masses; 7, , are the 
relaxation times; and el, g are the cyclotron 
frequencies. The plasma frequencies are 


51,27 [Ln, ,e*/m 12“ +xL)}, 

where n, , are the number of particles of the first 
and second kind, x is the electric susceptibility 
of the crystal in the absence of the free carriers, 
and L is the depolarization or sample shape fac- 
tor. The electric field is E= E,/ (1+ xL), where 
E, is the applied field (which is taken to be along 
a principal axis of the sample). It is assumed 


that both types of particles have the same charge e. 


For the interpretation of the present experiment 
we assume that the light- and heavy-hole masses 
are independent of energy and momentum. For 
the light holes this is an excellent approximation; 
the cyclotron mass of the heavy holes, however, 
has a momentum dependence which results in an 
anisotropy in the position of the cyclotron reso- 
nance of approximately 30%. We have used the 
appropriate cyclotron mass of the heavy holes 
for H, in the [110] direction in the above equations 
but have determined the ratio n,/n, = 20 from the 
density-of-states mass. 

In Fig. 1 the positions of the resonance lines, 
as calculated from Eqs. (1), are plotted as a func- 
tion of the total carrier concentration N. The ar- 
rows show the positions of the lines at low carrier 
concentration. The dashed lines indicate what 
would be the motion of each resonance if there 
were no coupling between the plasmas. In Fig. 1 
it is assumed that the microwave field is circu- 
larly polarized. For linear polarization the two 
directions of H, are equivalent, and Fig. 1 would 
have reflection symmetry in the horizontal axis. 
The thickness of the lines is a qualitative indica- 
tion of the relative peak intensity at a given car- 
rier concentration of the two resonances. It is 
thus seen that the absorption line which was in- 
itially the light-hole cyclotron resonance moves 
through zero field, increases in intensity, and 
eventually becomes the main magnetoplasma line; 
that which was initially the heavy-hole line de- 
creases in intensity, never passes through zero 
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FIG, 1. The calculated resonant magnetic field for 
the magnetoplasma resonances vs carrier concentra- 
tion. The thickness of the lines is approximately 
proportional to the relative intensity. The dashed 
lines indicate what would be the resonant fields if 
there were no coupling between the plasmas. The 
circles are experimental points, 


field, and asymptotically approaches a field close 
to the cyclotron resonance field for the light holes. 
It is clear that the absorption at each resonance 
cannot be ascribed either to light or to heavy holes 
but is due to the cooperative effect of both car- 
riers. The circles in Fig. 1 are experimental 
determinations of the resonance maxima positions. 
The discrepancy in the initial motion of the heavy- 
hole line is probably due to the assumption of iso- 
tropic mass. 

The weak resonance which is observed at low 
magnetic field for large carrier concentrations 
is of particular interest. In the present case, 
where wy9=1.58 Wp1> this resonance approaches 
4 the light-hole cyclotron resonance field for 
®p1,2>w. If W2>> Wp > w, the line would ap- 



























pear at the light-hole cyclotron resonance field. 
This resonance is an example of the phenomena 
discussed (for a somewhat different geometry) by 
Kittel’ in which it is possible to observe the cyclo- 
tron resonance of a minority group of carriers 

in the presence of a majority group because the 
latter group acts to short-circuit the plasma ef- 
fects of the first. This resonance is also similar 
to that discussed by Lax and Mavroides* for two- 
carrier systems in metals. 

Since the signal was not always a very small 
perturbation of the cavity, it was simpler to take 
the data with linearly rather than with circularly 
polarized microwaves, and this was done in spite 
of the fact that the two lines were at times super- 
imposed. In Figs. 2(a), (b), (c), and (d), the theo- 
retical absorption curves for linear polarization 
are plotted for several representative total car- 
rier densities N. The relaxation time for each 
group of carriers was chosen empirically to fit 
the experimental data.® The light solid line is the 
absorption due to the light holes, the dashed line 
that due to the heavy holes, and the bold solid line 
is the resultant. The arrows indicate the positions 
of the light- and heavy-hole cyclotron resonances 
at low carrier concentration. In Figs. 2(e), (f), 
(g), and (h) are experimental curves of the absorp- 
tion at approximately corresponding N. In Figs. 
2(a), (b), (e), and (£) the two lines are moving to- 
ward lower field with increasing carrier concen- 
tration. In Figs. 2(c) and (g) the strong line has 
passed through zero field and is moving to higher 
field as the carrier concentration increases. The 
weak line at about 0.8 koe is largely masked in the 
tail of the strong. At still higher concentrations 
[Figs. 2(d) and (h)] the main magnetoplasma line 
has moved off to the right and its tail is small 
enough that we can readily see the weak resonance 
line. The upper curve in Fig. 2(h) is the observed 
signal when the sample is placed in a regicn of 
uniform electric field in the microwave cavity. 
This is in approximate agreement with the theo- 
retically calculated curve. When the sample is 
placed close to a node in the electric field, the 
line is greatly enhanced with respect to the main 
magnetoplasma resonance, as shown by the lower 
curve. A simple physical argument indicates that 
at least part of this enhancement results from less 
cancellation of the currents due to the light and 
heavy holes when there is a gradient in the elec- 
tric field across the orbit of the carriers in the 
magnetic field. 

In many resonance experiments the sharpness 
and resolution of the lines can be greatly improved 
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FIG, 2. The upper set of curves shows the theoretical power absorption as a function of magnetic field for the 
carrier concentrations indicated. The dashed line is the absorption by the heavy holes, the light solid line is 
that by the light holes, and the bold line is the resultant. The bottom set of curves is the experimental absorption 
for the approximately corresponding number of carriers. The arrows indicate the position of the cyclotron reso- 
nances at low carrier densities. 


merely by making the observations at a higher ex- his assistance throughout the experiment. 
perimental frequency. This is not true in the pres- 





ent case. If we approximate the scattering rate 1G, Dresselhaus, A. F. Kip, and C. Kittel, Phys 
due to ionized impurities by aNT™ and the scat- , 100, 618 (1955). dit te a i , ‘ 
tering rate due to thermal phonons by 67**, where 2R, E, Michel and B. Rosenblum, Bull. Am. Phys. 
T is the temperature and a and bd are constants, Soc. 5, 375 (1960). 


and assume that these scatterings are independent, 5C, Kittel, Proceedings of the Conference on Radio 
and that neutral impurity scattering is negligible, and Microwave Spectroscopy, Duke University, 1957 


then there is a maximum possible w7 which can agree ss aici ducati 
be obtained for a given fractional shift of the reso- . and J. 5. Mavroides, in Solid ate Eaves, 
g . es 2 edited by F, Seitz and D. Turnbull (Academic Press, 





nance. This maximum is obtained by properly Inc., New York, 1960), Vol. 11. 

choosing the doping of the sample. The w7 in the It may be noted that the relaxation times for im- 

present work is not far from that maximum. purity scattering of the light and heavy holes are not 
We would like to thank Mr. Harold Hanson for equal while those for phonon scattering are. 
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ULTRASONIC AMPLIFICATION IN CdS 


A. R. Hutson and J. H. McFee 
Bell Telephone Laboratories, Murray Hill, New Jersey 


and 


D. L. White 
Bell Telephone Laboratories, Whippany, New Jersey 
(Received August 28, 1961) 


We have observed substantial amplification of 
ultrasonic waves in photoconductive CdS, pro- 
duced by applying a dc electric field in the direc- 
tion of wave propagation. Results are reported 
showing gains of 18 db at 15 Mc/sec and 38 db 
at 45 Mc/sec in a 7-mm length of CdS. 

Light-sensitive ultrasonic attenuation has been 
observed in CdS.*~* It was pointed out® that these 
effects result from the relaxation interaction of 
mobile charge carriers with the strong longitu- 
dinal electric fields of piezoelectric origin ac- 
companying certain acoustic waves in CdS. In 
the course of developing a theory for the propa- 
gation of acoustic waves in piezoelectric semi- 
conductors,® one of us (D. L. White) pointed out 
the possibility of achieving acoustic gain by ap- 
plying a de electric field which causes the inter- 
acting charge carriers to drift in the direction 
of wave propagation faster than the sound veloc- 
ity. 

Figure 1 shows the experimental arrangement. 
A 1-psec pulse of rf is applied to the top Y -cut 
quartz transducer producing a shear wave which 
propagates down through the first buffer, CdS 
crystal, second buffer, and is converted by the 
bottom transducer to an rf output signal. The 
output is amplified, detected, and displayed on 
an oscilloscope whose time base is triggered at 
the initial rf input pulse. (The fused silica buf- 
fers provide time delay and electrical insulation, 
The CdS crystal is oriented so that the shear 
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FIG, 1. Experimental arrangement. 





wave propagates at right angles to the hexagonal 
axis with material displacement along the hexag- 
onal axis. In the absence of electrical conductiv- 
ity, a longitudinal piezoelectric field proportional 
to strain (-e,,S/€) accompanies the wave in the 
CdS. Mobile charge carriers are produced with 
relative uniformity throughout the sample by illu- 
minating it with the weakly absorbed 5770/5790 A 
lines of mercury. The parallel polished faces of 
the CdS, which are bonded to the buffers, were 
prepared as electrical contacts by diffusing indium 
into a very thin surface region and then evaporat- 
ing a metallic indium film onto the surface. A dc 
voltage pulse can be applied between these con- 
tacts to provide the carrier drift field. The volt- 
age, duration, and time delay (with respect to 
the initial rf pulse) of the drift-field pulse are 
separately controllable. At the start of this ex- 
periment, the crystal appeared to be insulating 
in the dark with an acoustic attenuation of less 
than 7 db at 45 Mc/sec and less than 2 db at 15 
Mc/sec.* 

We have measured the change in the output sig- 
nal with respect to that obtained in the dark with 
no drift field, with the following results: 

1. No illumination—drift-field pulse has no ef- 
fect upon output signal. 

2. No drift field—increasing illumination yields 
increasing attenuation and conductivity. 

3. Sample illuminated—5-ysec drift-field pulse 
has no effect unless time delay is such that drift- 
field pulse overlaps 3.5-usec transit time of 1l-usec 
ultrasonic pulse in the sample. When overlap oc- 
curs, the output signal is either increased or de- 
creased, maximum effect corresponding to com- 
plete overlap. 

4. Sample illuminated, drift-field pulse com- 
pletely overlapping sound pulse in sample-varia- 
tion of the drift-field strength produces changes 
in the output signal as shown in Fig. 2. (For 
those parts of the curves showing large gain, 
input attenuation was added in order to remain 
within the linear region of the effect.) 

The striking features of these results are the 
appearance of negative attenuation (gain) and the 
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FIG, 2. Observed attenuation as a function of drift 
field. Zero db represents attenuation with sample un- 
illuminated. Curve A: 15 Mc/sec, w,/w=1.2. Curve 
B: 45 Mc/sec, w,/w=0.24. Curve C: 45 Mcfec, We /w 
=0.21. (w/w values are average values for the whole 
sample computed from measured CdS conductivities .) 


crossover at ~700 v/cm. The drift field at cross- 
over corresponds to a drift velocity, vg, equal to 
the appropriate CdS shear wave velocity (v =2 x10° 
cm/sec) for negative carriers of uy = 285 cm?/volt 
sec. This drift mobility is in good agreement 
with the Hall mobility of about 300 cm?/volt sec 

at 300°K reported by Piper and Halsted.” The 
change from acoustic loss to gain at a carrier 
drift velocity equal to the wave velocity corre- 
sponds to White’s prediction. 

A small signal theory of this attenuation behav- 
ior may be obtained by adding a constant drift 
field, Ez, to the equation for current density in 
reference 6. Basically, the theory amounts to 
obtaining an expression for the amplitude of the 
wave-periodic electric displacement, D, in terms 
of the wave-periodic E field from Poisson’s equa- 
tion and the current-density-space-charge con- 
tinuity equation. This relation between D and E 
is then used to eliminate D from the piezoelec- 
tric equation, D=eS+eE, so that an expression 
for E in terms of the strain, S, and the appropri- 
ate piezoelectric e constant is obtained, which 
may be substituted in the elastic equation, T=cS 
-eE. The result is a complex elastic stiffness 
constant, c’, which contains all of the electrical 
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effects. The acoustic absorption constant, a, is 
obtained from the imaginary part of (c’)"*. If 
one is dealing only with a single carrier species, 
we find 
oe (w/e) 
” a ao 2, 2 \2 |? 
v 2ce |1+ (w/v) (l+w /w py) 


(1) 





where (e?/ce) is the square of the appropriate 
electromechanical coupling constant; we =a/e, 
the conductivity relaxation frequency; y =1 -vq/v; 
and wp =v?/D, where D is the carrier diffusion 
constant. When the Einstein relation between mo- 
bility and diffusion is valid, (w?/w wp) =k’*A’, 
where k=27/) and A is the Debye length. Plots 
of the bracketed term in (1) are shown in Fig. 3. 
The effect of diffusion is to decrease the maxi- 
mum value of this term as w,/w is decreased. 
(If diffusion were neglected, the bracketed term 
would have maxima of +0.5 and the maximum 
values of +a would remain constant even in the 
limit of small w,/w.) The coefficient of the 
bracket term of (1) is ~ 150 db at 45 Mc/sec and 
= 50 db at 15 Mc/sec for our 7-mm sample." ** 
Thus, Figs. 2 and 3 have the same vertical scales 
and can be compared directly. 

Agreement between theory and experiment is 
semiquantitative for the data shown, which were 
taken at high light intensities. Closer agreement 
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FIG. 3. Plots of bracketed term of Eq. (1). 
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FIG. 4. Buildup of acoustic oscillations. Top trace: 
Drift-field pulse applied to CdS. Vertical scale, 715 
volt/em per large division. Lower trace: Signal from 
output transducer (transducer and receiver tuned to 15 
Mc/sec). Horizontal scale for both traces, 20 usec 
per large division. 


probably is not to be expected for the following 
reasons: Inhomogeneities in the photoconductiv- 
ity of the sample and a 7% ripple on the drift- 
field pulse probably caused blurring of the steep 
parts of the gain vs drift-field curves. Uncer- 
tainty may exist in the measured piezoelectric 
constant. There may be some photoconductivity 
due to holes. The Einstein relation (relating wp 
to 1) may fail because the electrons depart from 
thermal equilibrium under action of the drift 
field. 

To demonstrate unequivocally the existence 
of net acoustic gain in our sample, we have length- 
ened the applied drift pulse, applying no input 





signal. There is sufficient gain in the downward 
direction to offset the attenuation in the upward 
direction as well as the losses incurred on reflec- 
tion at the sample-buffer bonds, with the result 
that the sample breaks into oscillation. The lower 
trace in Fig. 4 shows the rapid buildup of acoustic 
energy from noise. A steady state is reached 
when the available carrier concentration can no 
longer support space-charge bunching sufficient 
to provide surplus gain at the high-strain ampli- 
tude. 

We are grateful to T. B. Bateman and H. J. 
McSkimin for advice and assistance in ultrasonic 
techniques and to D. E. Collins, L. J. Heilos, and 
W. Pleibel for sample preparation. 
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SPLITTING OF ELECTRON SPIN RESONANCE LINES BY AN APPLIED ELECTRIC FIELD 


G. W. Ludwig and H. H. Woodbury 


General Electric Research Laboratory, Schenectady, New York 
(Received August 17, 1961) 


Nuclei or paramagnetic ions in many solids oc- 
cupy sites which lack inversion symmetry. Bloem- 
bergen has recently called attention to the possi- 
bility of observing shifts, proportional to the ap- 
plied electric field &, in the energy levels of such 
systems.’ Consistent with Bloembergen’s ideas, 
Kushida and Saiki* and Armstrong et al.* have de- 
tected a broadening of certain nuclear quadrupole 
resonance lines. In this Letter we report the first 
observation of changes, linear in &, in the spin 
Hamiltonian of paramagnetic ions: The g tensor 
has been modified and its symmetry lowered, and 
a zero-field splitting of spin levels has been intro- 
duced. In the following Letter* a mechanism for 
such effects is proposed. 

The measurements were made on transition 
metal ions occupying interstitial sites in the sili- 
con lattice.® It was anticipated that effects of an 
electric field might be observable since the fol- 
lowing criteria are met: Interstitial sites lack 
inversion symmetry; the ¢, functions of the transi- 
tion metal ions are so hybridized that there are 
nonvanishing matrix elements of the electric field; 
and the separation between the ground and excited 
states of the ions are small. It is noted that some 
ions have orbital degeneracy which is partially 
lifted by spin-orbit interaction, and spectra of 
several ions are sensitive to uniaxial stress. 

Samples were studied with His using a conven- 
tional TE,,, mode cylindrical cavity. An &~12 
kv/cm could be applied along the long dimension 
of the sample. 

The ground state of interstitial neutral iron in 
silicon has an electron spin of one. So long as 
the tetrahedral symmetry is undisturbed, the 
three spin levels in a magnetic field are equally 
spaced, and the M =1 to 0 and 0 to -1 transitions 
are equal in energy. For &#0 the single reso- 
nance line splits, proportional to &, into two 
equally intense components, as shown in Fig. 1. 

It is believed that ions occupying one type of inter- 
stitial site display two fine structure components 
which differ in intensity by the Boltzmann factor 
(~1.17 under our experimental conditions). How- 
ever, there are two types of interstitial sites, 
differing by the operation of inversion, which 

give complementary spectra. The position of the 
double-quantum transitions® is not sensitive to 6, 
indicating that one is observing a fine structure 
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splitting rather thana g shift. For & ina cubic 
direction, the splitting varies sinusoidally for 
H.ié, as shown in Fig. 2. 

The ground state of interstitial Fe* in silicon 
can be described by an effective total angular mo- 
mentum J= 4; there are two energy levels and 
one resonance transition.* For & sufficiently 
large one observes two equally intense resonance 
lines, which show an angular variation similar 
to that described for Fe®. Here, however, the 
appearance of two lines must represent an effect 
on the g tensor; the shift is of opposite sign for 
the two nonequivalent sites. 


kv/cm 





10 
kv/cm 





t 


FIG, 1. Splitting of the spectrum of Fe’ in silicon by 
an electric field applied in the [001] direction, for Hi in 
a [110] direction. The temperature is 4, 2°K. 
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The spin Hamiltonian describing the results 
for Fe° and Fe* with g parallel to the [001] direc- 
tion is 
K= g8J-H+ &[68(J,H, -JgH,) +E, -J,")). (1) 


Here the subscripts 1 and 2 denote the [110] and 
the [110] directions, respectively, and the plus 
and minus signs correspond to the two different 
sites. To first order in € the M to M-1 transi- 
tions are 


hv =g8H+ &[bBH +3E(M -}4) |sin*6 cos2¢9, (2) 


where @ and g are polar angles describing the 
direction of H, and 6=7/2, g=0 is the [110] direc- 
tion. We find b~ 0 and |E| =(4.7+ 0.2) x107~°/kv 
for Fe°®; |b) =(7+1) x10 cm/kv for Fe*. Equa- 
tions (1) and (2) are easily generalized to arbi- 
trary direction of &.* 

Effects of electric fields have also been ob- 
served for other ions in silicon. 

One of us (G. W. L.) wishes to acknowledge his 
indebtedness to Professor Bloembergen, who sug- 


30 60 
DIRECTION OF MAGNETIC FIELD (¢) 


gested the possibility of observing electric field 
effects in silicon in an informal discussion. We 
are grateful to F. S. Ham for discussions which 
contributed to this work. C. R. Trzaskos assisted 
in the measurements. 
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LINEAR EFFECT OF APPLIED ELECTRIC FIELD IN ELECTRON SPIN RESONANCE 


Frank S. Ham 
General Electric Research Laboratory, Schenectady, New York 
(Received August 17, 1961) 


In the preceding Letter,’ Ludwig and Woodbury 
report observation of linear effects of an applied 
electric field on the electron spin resonance 
spectra of two charge states of interstitial iron 
in silicon, in accord with Bloembergen’s predic- 
tion” that such effects should be in general observ- 
able for ions at crystal sites lacking inversion 
symmetry. The purpose of this Letter is to show 
explicitly how these effects arise for iron-series 
transition metal ions that form interstitial defects 
in silicon. Similar results pertain to other ions 
with d” or f” valence configurations in various 
crystals. 

A phenomenological model that successfully 
describes interstitial iron series ions in silicon 
has been given by Ludwig and Woodbury*: The 
ions occupy sites of tetrahedral symmetry; all 
valence electrons are in the d shell; the three- 
fold degenerate ¢, orbitals are lower in energy 
than the doubly degenerate e’s; Hund’s rule is 
obeyed. However, the ¢, and e wave functions 
of the interstitial are very different from the 3d 
functions of the free ion,* although they have the 
same behavior under symmetry operations of 
the tetrahedron. In particular, 3d¢t, functions 
mix with functions of odd parity on the metal ion, 
such as 4p, and form bonds with neighboring 
silicon atoms. Denoting by x, y, z coordinates 
with respect to the cube directions, we may rep- 
resent the /, functions schematically by 


v =a(yz)5 +B), +79, “Py +P -@ 4), (1) 


and Yy and wz by cyclic permutation of x, y, z. 
Here Y,, 2, Y;, and g, denote o orbitals on 
nearest neighbors in the [111], [111], [11I], and 
[111] directions, respectively, and for simplicity 
we ignore other possible bonds. 

The interaction Kg =e 6,2 of an electron with 
an electric field & 2 in the z direction then has 
only one nonzero matrix element between the ¢, 
functions: 


leg zl - : 2 
(y é, v, ps. (2) 
This defines a parameter p. We may for conven- 
ience represent the matrix elements between f, 


functions in terms of an operator £ which satis- 
fies the commutation rules for angular momentum, 
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with £=13: 


Ke=-E WL £ +2, £). (3) 


y 
HK g lifts the orbital degeneracy of ¢,, energy 
levels being at 0 and+p&,. Thus an electron 

in a ¢t, state in tetrahedral symmetry responds 

to a field in a cube direction as if the state pos- 
sessed an electric dipole moment of strength p. 
From (1) and (2), p is nonzero if the ¢, states 
are mixtures of functions of opposite parity on 
the central ion (a@ and B# 0) or if o bonding oc- 
curs with the nearest neighbors (y #0). 

For Fe® or Mn” interstitial in Si, the ground 
state arises from the configuration ¢,°e? and has 
S=1 and no orbital degeneracy.* At an energy A 
higher, there is an orbital triplet of the form 
t,*e* which is split by an electric field. Between 
these states and the ground state there are ma- 
trix elements of the spin-orbit interaction r(L-5), 
the orbital part of the Zeeman interaction B(L- 8), 
and the interaction P(LI) between the nuclear 
moment and the orbital moment of the electrons. 
Accordingly, from third-order perturbation the- 
ory we find for the spin Hamiltonian of the ground 
state® 


K=gBH-S+Al-S -78 yHeT 
+ , (2.?pa?/A?) ee “8 
+ 6 pagan" (EOF, +8 8) 
+6 (4apPaa”’ /?) St + St) 
+ 6 (2pPta’?/a*\(I 1 D1) 
+& By (4pPa ra" /)(B By) tH +1 i )- (4) 


Terms in (4) arising from &,. and &, may be ob- 
tained by permutation of indices. Here a, a’, 
and @’’ are factors arising because matrix ele- 
ments of A(L- 8), 8 (L+H), and P(L-f) are changed 
from their free-atom values. The isotropic shift 
in g provides the relation 


(8raa’/A) = -(g - 2.0023). 


For gBH > (\*a?/a*)pé and arbitrary orienta- 
tion of 6 and H, the energy levels to first order 
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in & are given by (with omission of nuclear terms 
E(M) =gBHM + (6A*pa?/A?) &M?o 
+(8rpaa’/A?)&§ BHM 6, (5) 
where 
o= ee +7 AS, ts a (6) 
Ny, Ny, nz and Mx, My, Mz being the direction 
cosines of H and &, respectively. 
For Fe* or Mn’ the ground state is ¢,5e? withS 
= 3 and a threefold orbital degeneracy describable 
by an effective orbital angular momentum opera- 
tor £ with £=1. With J=£+S, the spin-orbit 
interaction separates the levels with J=}, 3, and 
3, placing J=4 lowest. The effect of &, is rep- 
resented by an operator of the same form as (3), 
and this, the Zeeman interaction 


> =g BRS H-S 
8p BH-£ +g HS, 
and the hyperfine interaction 
=A £.T ST 
iH, 'e 1+A 5 I 


have matrix elements between the J= 4 and 3 
states. Accordingly, from second order pertur- 
bation theory the spin Hamiltonian for the J=4 
states has the form 


affie « H.J4+(8A -2 mi 
Ki (3g. ig , BH J+(3A . §A ,)J-l 


‘ [2p(g. -g 2 3W) & 8 UA, +3, H) 
-[2p(A , -A,)/3W] 8, 1, +I1), 


where W =[E(%) -E(4)]. 

For interstitial Mn* and Cr° in silicon the 
ground state is t,*e? (S=2, £=1). The effect of 
&, is given by (3). This gives a direct Stark 
splitting of the ground level (J=1) as well as 
producing shifts in other resonance parameters. 

The silicon lattice has two tetrahedral inter- 
stitial sites identical in their environments ex- 
cept that one is inverted with respect to the other. 
The dipole moment p has opposite sign for these 
sites, and resonance spectra should thus show a 
symmetrical splitting of lines to first order in 6 , 
as is observed. 


We have so far ignored the fact that the para- 
magnetic ion is displaced by the external field 
from its position of tetrahedral symmetry. This 
displacement gives rise to a perturbation in the 
potential seen by the electrons in addition to that 
produced directly by the applied field. However, 
we may show from symmetry that to first order 
in & this further perturbation can affect only 
those matrix elements already affected directly 
by the applied field. Thus the parameter p will 
be used to represent the net effect of both per- 
turbations. 

From the results of Ludwig and Woodbury and 
Eqs. (4) and (7) we can evaluate the effective one- 
electron dipole strength p for interstitial Fe* and 
Fe® in silicon, using estimates of the other param- 
eters in these formulas obtained from other stud- 
ies.* We obtain p/e =0.3 x10 cm for both Fe® 
and Fe+. These values are, however, each un- 
certain to perhaps a factor of 3 because of uncer- 
tainties in making correction for local fields and 
in estimating the other parameters. 

If p is of the same magnitude for Mn and Cr in 
silicon as for Fe, the ground states of Mn* and 
Cr° should suffer a linear Stark splitting of ~0.05 
cm™ with §=10* v/cm, in zero magnetic field. 

The author is indebted to G. W. Ludwig and 
H. H. Woodbury for helpful discussions. 
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DE HAAS-VAN ALPHEN EFFECT IN RHENIUM, NIOBIUM, AND TANTALUM* 


A. C. Thorsen and T. G. Berlincourt 


Atomics International, Division of North American Aviation, Inc., Canoga Park, California 
(Received July 24, 1961; revised manuscript received August 21, 1961) 


The electronic structures of many low-melting- 
point metals have been extensively investigated 
and characterized by de Haas—van Alphen (dHvA), 
galvanomagnetic, cyclotron resonance, anomalous 
skin-effect, and magnetoacoustic absorption tech- 
niques.’ In contrast, despite considerable inter- 
est in the electronic character of transition met- 
als,” there have been very few applications of 
such techniques to these elements*’* due primarily 
to the difficulties of obtaining single crystals of 
the requisite high purity. In this paper we report 
on the first observations of dHvA oscillations in 
the transition metals Re, Nb, and Ta. As will be 
evident from the data, the ranges of the observed 
dHvA parameters are similar to those observed 
for polyvalent low-melting-point metals. 

The single crystals of Re and Nb used in this 
investigation were grown, respectively, by Nad- 
ler and Buehler using electron bombardment 
techniques. The Ta sample was prepared by 
Budnick using a strain-anneal technique.’ Chase 
Brass and Copper Company material was used 
for the Re crystal which had a room-temperature 
to liquid-helium-temperature resistivity ratio of 
540. The corresponding ratios for the Nb and Ta 
crystals were 56 and = 600, respectively. (The 


liquid helium temperature resistivities were 
measured at 4.2°K in magnetic fields just suf- 
ficient to quench all superconductivity.) The Re 
and Nb samples were approximately 1 mm in di- 
ameter by 9 mm in length and were obtained from 
single-crystal rods by a combination of etching 
and spark cutting techniques. The Ta sample 
was cut from a 0.51-mm diameter wire and was 
composed of several small single crystals whose 
orientations varied slightly along the length of 
the sample. The approximate crystallographic 
orientations for all three samples are indicated 
in Table I. The pulsed-field apparatus used in 
this study has been described briefly.**” 

A complex electronic structure might well be 
anticipated in Re. If one assumes seven valence 
electrons per atom, the free-electron Fermi 
sphere completely encloses the first two Bril- 
louin zones for the hexagonal lattice. In such a 
situation, there may exist in higher zones a num- 
ber of electron and hole surfaces capable of gen- 
erating dHvA oscillations. The results are sug- 
gestive of such a complex structure, for no fewer 
than four oscillating susceptibility terms were ob- 
served in Re with the field approximately parallel 
to the [1010] direction. Three of these terms can 


Table I. Periods, corresponding Fermi surface extremal cross-sectional areas, and effective masses for 
de Haas— van Alphen oscillations in Re, Nb, and Ta. (Estimated accuracy of periods is +5 %.) 








Magnetic field 


P A 
(10'5 em=*) 





Crystal orientation (1078 gauss~') m*/my 
Re [1010] 1.65 5.79 1.9 +0.7 
6.3 1.52 0.5 +0.05 
7.1 1.34 0.53 +0.05 

124 0.077 <0.1 

Nb [110] 1.16 8.23 1.0 +0.1 
1.67 5.72 1.0 +0.1 

Ta® 12° from [110] 2.02 4.73 

in (001) plane 2.16 4.42 0.8 +0.1 


————- 





a 
Only one effective-mass entry appears after Ta because the required data were taken during a second run when 


only one term was observable. 
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FIG. 1. de Haas—van Alphen effect in Re at 1.16°K. FIG, 2. Long-period de Haas—van Alphen oscilla- 
The oscillating trace shows the output from a pickup tions in Re at 1.16°K. The diagonal trace shows the 
coil containing the sample. The curved trace shows field increasing (from right to left) from 13.5 to 42.1 
the field increasing from 87.7 to 94.4 kilogauss during kilogauss during a sweep time of about 2.0 milliseconds. 


a portion of the oscillating trace. The total sweep time 
is two milliseconds (time increasing from right to left). 
Three oscillating terms are visible. ing terms exist but were masked by collision 
broadening of the Landau levels. The relatively 
low resistivity ratio suggests that such broaden- 
be recognized in Fig. 1, while the fourth and much _ ing might be quite pronounced. Little can be said 
longer period term, observable at lower magnetic as yet regarding the origin of the dHvA oscilla- 
fields, is shown in Fig. 2. The periods P deduced tions in Nb except that those we observe are prob- 
from such data are given in Table I along with the ably not due to electrons or holes in the first two 
corresponding Fermi surface extremal cross-sec- bands. The free-electron sphere completely en- 
tional areas A =417e/chP. The effective-mass closes the first Brillouin zone for a pentavalent 
values m* listed in Table I were calculated from bec metal, and on the basis of a nearly free elec- 
the temperature variations of the amplitudes of 
the oscillations in the usual manner.* It should 
be remarked, however, that for the longest peri- 
od oscillations the corresponding m* is so small 
that the exponential approximation in the temper- 
ature dependence of the amplitude was inapplica- 
ble and a comparison with the hyperbolic sine 
function was required. Even so, the temperature 
dependence of the amplitude was such that only 
an upper limit could be assigned to the effective 
mass. 

The dHvA parameters in Re have values typical 
of those for nontransition metals. In fact, the 
longest period oscillations observed are compara- 
ble with those which occur in semimetals and 
should be easily observable with greater precision 
in modest magnetic fields by means of torsion 
balance techniques. 

Two oscillating terms were observed in Nb 








FIG. 3. de Haas—van Alphen effect in Nb at 1.14°K. 
The beating effect is due to two oscillating terms of 






with the magnetic field approximately parallel different period and amplitude. The diagonal trace 
to a [110] direction. These are pictured in Fig. shows the field decreasing (from right to left) from 
3, and the corresponding parameters are listed 150.0 to 138.2 kilogauss during a sweep time of 1.6 






in Table I. It is probable that additional oscillat- milliseconds. 
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FIG. 4. de Haas—van Alphen effect in Ta at 1.13°K. 
The beating effect due to two oscillating terms can be 
seen as the field decreases (from right to left) from 
141.6 to 134.9 kilogauss during a sweep time of 1.0 
millisecond. 


tron model®’® the holes in the second band would 
be expected to give rise to periods shorter than 
those observed. 

Two oscillating terms were also observed in Ta 
(see Fig. 4), and the corresponding dHvA param- 
eters are listed in Table I. Our measurements 
indicate that the shorter period term is quite sen- 
sitive to orientation, for it nearly disappeared 
when the sample was removed and remounted at 
a slightly different angle. It is interesting that Ta 
and Nb have the same valence and crystal struc- 
ture, and have lattice parameters differing by 
less than one percent. For such a situation the 


nearly free electron model predicts almost iden- 
tical periods. The observed differences (less 
than 50%) could be attributed to the difference 
in crystallographic orientation. 

The authors wish to thank E. Buehler and J. E, 
Kunzler of Bell Telephone Laboratories for sup- 
plying the sample of Nb, J. I. Budnick of IBM 
Watson Research Laboratory for the Ta sample, 
and H. Nadler of this laboratory for the Re sam- 
ple. We would also like to acknowledge many 
valuable discussions with R. E. Behringer. 
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Paramagnetic resonance techniques have been 
used to measure exchange interactions between 
nearest neighbor (nn) Mn** ions in mixed crys- 
tals of Mn, ZnF, with Mn:Zn~1:50. The two ions 
are separated by the short edge (c ~ 3.2 A) of the 
body -centered tetragonal unit cell. The general 
procedure has followed that used for previous in- 





S'= SJ = §. 
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R. W. H. Stevenson 
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vestigations'»* of exchange-coupled pairs of spins, 


FERROMAGNETIC EXCHANGE INTERACTION BETWEEN Mn”* IONS IN Mn, ZnF, 


M. R. Brown, B. A. Coles, and J. Owen* 
Clarendon Laboratory, Oxford, England 


The results can be described by the following 
spin Hamiltonians. For one ion of the pair, omit- 
ting interactions with the other ion, 

sc’ = got 3' + [8,')” - 45's" +1)] 
i,2 i,2 ( 
e 1) 
tE IY - 61, 


where g=2, S'=$, D, = -0.0135+ 0.004 cm™, E¢ 
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=0+0.006 cm™. There is a similar Hamiltonian 
x for SJ. For the interacting pair the Hamil- 
tonian is 


sd aad +30! + 8'.84 «0 [5,'s 7-384, @) 


where J/k = -0.4+ 0.2°K and D, = -0.051 + 0.009 
em™, The z axis is parallel toc, that is, to the 
line joining the ions. The isotropic exchange 
term JS’. §/ gives rise to states of total spin 
(§=8'+8/) with S=0, 1, 2, 3, 4, and 5, and with 
relative energies, 0, J, 3J, 6J, 10J, and 15J, 
respectively. Since J is negative for the present 
case, S=5 is lowest and S=4 is 5|J| higher in 
energy, and most of the measurements have been 
made on transitions within these two lowest total 
spin states. 

The values of the parameters D- and De, which 
determine the position and anisotropy of the lines 
in the spectrum, have been confirmed by fitting 
more than 20 transitions in the wavelength range 
0.85 to 3.3 cm. The crystal-field contributions, 
D,= -0.0135 cm™, Ec=0, are of the same order 
as those found for isolated Mn** ions in very di- 
lute crystals which are* D, = -0.0186 cm™, E, 
=-0.0041 cm™. The other contribution Dg ap- 
pears to be consistent with the expected dipole- 
dipole interaction -g*/r;;*, since the value can 
be fitted by choosing A he $01 A which lies be- 
tween the expected* separations 3.13 A for ZnF, 
and 3.31 A for MnF,. 

The magnitude of the isotropic exchange J has 
been estimated by using intensity measurements 
to find the separation between different energy 
levels given by Eq. (2). This has involved the 
determination of the intensities of a number of 
lines in the spectrum both relative to each other 
and relative to a comparison specimen in the 
temperature range 1.5 to 20°K. All the meas- 
urements agree in giving the exchange ferro- 
magnetic in sign and close to the value Jpy/k 
=-0.4°K. Preliminary measurements have also 





been made on next nearest neighbor (nnn) Mn- Mn 
pairs which give the approximate result Jpnp/k 
= 4°K (antiferromagnetic). 

This preliminary result for next nearest neigh- 
bors is in reasonable agreement with a molecular 
field interpretation of the bulk properties of the 
antiferromagnetic MnF, as was expected. The 
ferromagnetic sign of the nearest neighbor ex- 
change, on the other hand, was not expected. If, 
as seems reasonable, a similar ferromagnetic 
interaction is present in MnF,, then in order to 
account for the susceptibility properties it may 
be necessary to assume appreciable antiferro- 
magnetic exchange between further neighbors on 
the same magnetic sublattice, e.g., third, fourth, 
sixth, or seventh nearest neighbors. 

Finally, it is of interest to note that the nearest 
and next nearest neighbor pair structure in MnF, 
is in some ways similar to that in MnO but the 
ratio Jpnn/Jnn is surprisingly different, being 
about -10 for MnF, and +1 for MnO.’ This is not 
understood, but a possible explanation may be that 
there is direct ferromagnetic exchange of magni- 
tude a few degrees Kelvin between nearest neigh- 
bors in both salts, and that this is the dominant 
interaction in MnF, where the antiferromagnetic 
superexchange interactions are smaller. 
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ish Board of Admiralty for supporting this work, 
and we are grateful to Dr. W. Marshall for help- 
ful discussions. 
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POLARIZATION IN NEUTRON-PROTON SCATTERING BELOW 100 Mev 


P, H. Bowen, G. C. Cox, G. B. Huxtable, A. Langsford, J. P. Scanlon, and J. J. Thresher* 
Atomic Energy Research Establishment, Harwell, Berkshire, England 
(Received August 17, 1961) 


A systematic investigation of neutron-proton 
scattering in the energy range 20 to 120 Mev is 
being made by this group with the object of pro- 
viding sufficient information for a phase-shift 
analysis of the T=0 state of the two-nucleon sys- 
tem. Since no polarization data on the n-p system 
have hitherto been available below 77 Mev, our 
first results appear to be of sufficient interest to 
publish at this stage. 

A polarized neutron beam has been obtained 
from the neutron time-of-flight spectrometer of 
the Harwell cyclotron’’?; the beam polarization 
as a function of energy is given in Table I. Fig- 
ure 1 is a schematic diagram of the apparatus, 

A liquid hydrogen target was placed at the end of 
the neutron flight path, and neutrons scattered 

by the target were detected by two large organic- 
liquid scintillation counters, placed symmetrically 
above and below the beam. In order to reject 
events in which a recoil proton from the hydrogen 
target was detected, anticoincidence counters 


Table I, Polarization of the neutron beam as a func- 
tion of energy. 





Neutron energy (Mev) Neutron beam polarization 





22.5 0.135 +0.029 
30 0.139 +0.014 
40 0.130 40.013 
50 0.187 +0.014 
60 0.262 +0.016 
70 0.312 +0.018 
80 0.338 +0. 022 
90 0.367 +0.029 
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were placed in front of the neutron counters. The 
polarizations were determined by measuring asym- 
metries in neutron scattering and, to reduce false 
asymmetries, a solenoid was used to rotate the 
polarization vector of the neutron beam alternate- 
ly clockwise and anticlockwise through approxi- 
mately 90°. The use of a symmetrical two-count- 
er system made it possible to analyze the data in 
such a way that the calculated scattering asym- 
metry was independent of the neutron beam inten- 
sity and the efficiencies of the counters. A beam 
monitor was necessary only to enable corrections 
for backgrounds to be made. These were deter- 
mined in subsidiary measurements and were due 
to rescattering of neutrons from the walls of the 
experimental area after they had been scattered 
by the hydrogen target and to scattering of neu- 
trons from the thin metal windows of the target. 
At neutron energies above 80 Mev, the number 
of detected neutrons varied rapidly with time of 
flight, and a small uncertainty in the alignment of 
the time spectra from the two neutron counters 
could have produced a large uncertainty in the 
measured asymmetry; for this reason it was 
necessary to reject most of the data at neutron 
energies above 85 Mev. Recent improvements in 
the electronic apparatus have virtually eliminated 
this source of error for future measurements. 
The asymmetries were corrected for the effects 
of multiple scattering of neutrons in the hydrogen 
and its metal container. From preliminary esti- 
mates a correction factor which increases the 
asymmetry by about 13 % of its value at 22.5 Mev, 
falling to 4.4% at 90 Mev, has been applied. The 
possible systematic error in making this correc- 
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FIG, 1. Experimental arrange- 
ment—schematic diagram. 













tio: 
tud 
the 


err 
eac 
fro) 
side 


ban 
eac! 
wer 
banc 
spec 





61 


‘he 
3ym- 
lse 


ite- 


ated 
ects 
gen 


sti- 


fev, 
The 


nge- 








VoLUME 7, NUMBER 6 


PHYSICAL REVIEW LETTERS SEPTEMBER 15, 1961 









































4 
+ 
~ = un 22.5Mev | pal gh PF a 
PO ° 
7 
a 7 
4r 
a ; P.2+ t t 1 
] ty t+ ,  30Mev “Tl +t 70 Mev 
Po $ o 
-.2 4 
FIG, 2. Polarization in Sr } 7 
neutron-proton scattering. at : At 
edictions of t P 
The curves are predic ons 0 2 40 Mev - 2b 8O Mev 
n-p polarization at 18, 40, and P e ps 
77 Mev, made on the basis of | t 
the Signell-Marshak potential ~2 4 
model. 6 
4 . 4 | + 
P ia + 7 P > 
“Ty ytt't + SOMev ‘ 90 Mev 
” 2P ak 6O°BOPIOO® 140° 180° 20°40°60°BOIOO” 140 180° 
a -—— 8m, 


tion, which should not be larger than the magni- 
tude of the correction, has been combined with 

the somewhat larger statistical error in the meas- 
urement of beam polarization to give a “scale 
error” at each energy. The individual points on 
each curve have statistical errors which arise 
from the asymmetry measurements and are con- 
siderably larger than the scale errors. 

The results have been grouped into 8 energy 
bands and are shown in Fig. 2 and Table I. At 
each energy, the over-all energy resolutions 
were obtained by combining the width of the energy 
band with the intrinsic energy resolution of the 


agreement with those obtained at 95 Mev* and 
77 Mev ® with considerably worse energy resolu- 
tion. 

Some predictions for n-p polarization, based 
on the Signell and Marshak potential model, are 
shown in Fig. 2 at 77 Mev® and at 18 and 40 Mev.°® 

The results presented here have enabled Phillips’ 
to place some restrictions on proton-proton scat- 
tering phase shifts at low energies, assuming 
charge independence of nuclear forces. 

A recent phase-shift analysis for the neutron- 
proton system,® which did not include the results 
presented here, has provided a number of possible 








spectrometer. The results are in reasonable sets of phase shifts, all of which give a good fit 
Table II. Neutron-proton polarization as a function of energy and angle. 
Total 

Nominal energy Scale 

energy spread error Neutron-proton polarization 

(Mev) (Mev) (%) 20° 30° 40° 50° 60° 80° 
22.5 5 25 0.03+0.09 -0.0320.12 -0.02+0.10 0.20+0.08 0.10+0.09 -0.12+0.09 
30 10 16 0.13 +0.06 0.05+0.07 -0.04240.06 0.11+0.04 0.08 +0.06 0.05 +0. 04 
40 10 12 -0.07 +0.07 0.14 +0. 08 0.30+0.08 0.17+0.05 0.18 +0.07 0.15 +0. 05 
50 10 10 0.09 +0.05 0.16 +0. 06 0.18+0.06 0.21+0.04 0.16 +0.06 0.16 +0. 05 
60 11 8 0.14+0.04 0.21 +0. 05 0.09+0.05 0.23+20.03 0.23 +0.06 0.15 +0. 04 
70 12 8 0.14+20.05 0.14 +0. 06 0.1420.06 0.22+0.04 0.36 +0.08 0.28 +0. 05 
80 13 8 0.28 +0.06 0.27 +0. 08 0.2620.08 0.3720.05 0.59+0.12 0.34 +0. 08 
90 14 9 0.27 +0.11 0.43 +0.08 


a 
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to the present data.* Further measurements are, 
however, being made to reduce the statistical 
errors of the present data, and to extend the re- 
sults to 160° c.m. and to 100 Mev. 





*Now at the National Institute for Research in Nuclear 
Science, Harwell, Berkshire, England. 
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EVIDENCE FOR CHARGE INDEPENDENCE IN MEDIUM WEIGHT NUCLEI" 
J. D. Anderson and C. Wong 


Lawrence Radiation Laboratory, University of California, Livermore, California 
(Received July 26, 1961) 


The importance of isotopic spin considerations 
(AT=0) in “mirror nuclei” (p,”) reactions has 
been pointed out by Bloom et al.’ If one goes to 
nonmirror nuclei such as V™, Cr™, the (p,) re- 
action is assumed to go as follows: The incoming 
proton reacts with an f,, neutron, exchanges its 
charge, and is emitted as a neutron. In the initial 
state there are eight f,,. neutrons available and 
three f,. protons. Since, by definition of isotopic 
spin state (charge independence), all the nuclear 
interactions within the initial and final nucleus 
are the same, the Q for the (p,m) reaction be- 
tween V™ (T= §) and its analog state in Cr™ is 
the Coulomb energy difference. It is the purpose 
of this Letter to point out that direct-reaction 
neutrons from the (p,m) reaction on medium 
weight nuclei do indeed leave the residual nucleus 
in a state which is the analog of the ground state 
of the bombarded nucleus. 

It was pointed out previously’ that the energy 
and angular distribution of the continuum neutrons 
from the V™(p,n)Cr™ reaction are in agreement 
with predictions of the statistical model of the 
compound nucleus for proton bombarding energies 
up to 8 Mev. Using time-of-flight techniques*»* 
to measure the neutron spectra, direct-reaction 
neutrons are observed for proton bombarding en- 
ergies above 10 Mev. These neutrons do not come 
from groups leaving Cr™ in low-lying states as 
might be expected from a simple Born approxi- 
mation calculation in which one considers only 
the radial overlap of the wave functions describ- 
ing the initial and final states. These neutrons 
are observed to come from a level (or levels) in 
Cr™ at 6.5 Mev. 


250 


Neutron spectra at 67; = 23° have been measured 
for proton energies between 9 and 13 Mev using a 
self-supporting 8-mg/cm? vanadium foil and a 
10-meter flight path. The neutron spectra for 
three incident proton energies are shown plotted 
in Fig. 1. The compound statistical model pre- 
dicts that if a sufficient number of nuclear levels 
are involved in both the compound nucleus and 
the residual nucleus, the energy distribution of 
neutrons emitted is given by 


P(E)dE =KEw(E )o E)aE, 


excit 
where w(E ox, jt) is the level density of the resid- 
ual nucleus at its excitation energy, which is 
determined by the incident proton energy, the Q 
value of the reaction, and the emitted neutron en- 
ergy; 0,(E) is the reaction cross section for the 
inverse reaction between the excited residual 
nucleus and a neutron of energy E; and the con- 
stant K is a function of the incident proton ener- 
gy. Assuming that o,(£) is a slowly varying func- 
tion of neutron energy; the energy dependence of 
the level density of the residual nucleus is propor- 
tional to P(E)/E. The relative neutron spectra 
transformed to the center-of mass system and 
divided by the energy of the emitted neutron in 
the c.m. system are plotted as a function of the 
excitation of the residual nucleus in Fig. 2. From 
Fig. 2 it is clear that P(E)/E (allowing for the 
change in energy resolution as a function of neu- 
tron energy) is only a function of the excitation 
of the residual nucleus up to Eexcit= 6.5 Mev, i- 
dicating the assumption of compound nucleus for- 
mation in this region is valid. However, a fluctu- 
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FIG. 1. Neutron spec- 
tra from proton bombard- 
ment of V*!, 











ation in level density of a factor of 2 to 3 at Eexcit 
=6.5 Mev with a hundred levels being involved is 
incompatible with this level (or levels) being popu- 
lated by a compound nucleus. As one increases 

the proton bombarding energy the continuum neu- 
tron cross section (Eexcit= 6.5 Mev) decreases 





reo Ft 








b 





FIG, 2. The neutron spectra transformed to the 
center-of-mass system and divided by the energy of 
the emitted neutron are plotted as a function of the ex- 
citation of the residual nucleus Cr®!, The arrow de- 
hotes the position of the analog, state. 





from 17 mb/sr Mev at E,=10.8 Mev to 4.5 mb/sr 
Mev at E,=13.3 Mev, while the cross section for 
the “level” remains at 2 mb/sr. The decrease in 
cross section for the continuum neutrons is in ex- 
cellent agreement with our previous measurement 
of the Cr™ level density*»* and is to be expected on 
the basis of a compound nucleus while the direct- 
reaction cross section is expected to vary only 
slowly with energy which is borne out by the pres- 
ent experimental data. Thus it is concluded that 
the level at 6.5 Mev in Cr®™ is excited via a direct 
reaction. 

If one calculates the V™-Cr™ Coulomb energy 
difference using the results of Swamy and Green,® 
one obtains AE, = 8.2 Mev for R=1.25A™“ f. This 
is in excellent agreement with our measured value 
of 8.0+0.2 Mev. [This is the (p,m) ground state 
Q(-1.54 Mev) plus the residual excitation in Cr™ 
(6.5 Mev).] Our conclusion is that the level in Cr™ 
at 6.5-Mev excitation is the analog of the ground 
state of V™ and is therefore a T= § state. 

Additional evidence for analog states is found in 
the Fe®®(p,n)Co® and Co®*(p,n)Ni® neutron spec- 
tra. Our Q values for the (p,m) reaction to the 
analog states are again in excellent agreement 
(see Table I) with the Coulomb energy differences 
calculated from reference 5, 

Further experiments involving additional target 
nuclei and measurements of the angular distribu- 
tions of the direct-reaction group are currently 
under way. 

It is a pleasure to acknowledge the assistance of 
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Table I. @ values for the (p,”) reactions to the analog 
state. 





Q (measured) 4Eg (calculated) 





(Mev) (Mev) 
v54(p,n)Cr5t # 8.0+0.2 8.2 
Fe*(»,n)Co® * 8.9+0.2 9.0 
Co®*(p,n) Nid? * 9.1+0.2 9.2 





a 
Denotes isotopic spin analog state. 


J. W. McClure and B. D. Walker in obtaining the 


data and of S. D. Bloom in its interpretation. 
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PHOTODISINTEGRATION OF POLARIZED AND ALIGNED DEUTERONS* 


W. Zickendraht, D. J. Andrews, and M. L. RustgiT 
Yale University, New Haven, Connecticut 
(Received August 7, 1961) 


Information regarding dynamic and electromag- 
netic properties of the two-nucleon system obtain- 
able from the photodisintegration of the deuteron 
can be extended by employing aligned and polar- 
ized deuterons. In the present note a partial and 
preliminary survey of the possibilities is attempt- 
ed by confining oneself to unpolarized gamma rays 
and to differential cross sections without a consid- 
eration of the polarization of ejected nucleons. 
Since, as will be seen below, the effects of align- 
ment are quite pronounced, one may hope that 
the assumptions underlying the usual theories of 
the photodisintegration process can receive a 
more thorough test by comparison of theory and 
experiment even to the limited extent gone into 
here. The calculations reported on have been 
made by employing a previously described ar- 
rangement’ in terms of probability amplitudes 
Smi Corresponding to final states with combined 
proton-neutron spin projection m arising from 
initial states with spin orientation i. A modified 
Signell-Marshak potential, referred to as Poten- 
tial I in reference 1, was used at six photon ener- 
gies in approximation E in which all transitions 
of types E1, M1, and E2 are considered in the 
nonretarded approximation. 

The line of orientation axis determined by the 
direction of fields external to the deuteron is a 
convenient choice for the quantization axis. In 
thermal equilibrium one deals essentially with a 





statistical mixture of magnetic substates with in 
general unequal population probabilities and with- 
out correlation. The calculations reported on are 
concerned with this somewhat idealized situation. 
Denoting the relative population numbers of the 
magnetic substates » by w, with the normaliza- 
tion 


W,+W,o+w.,=1, (1) 


the polarization parameter P, and the alignment 
parameter P, in the notation of Blin-Stoyle and 
Grace’ are, respectively, 


P, =W,-W.4, (2) 
P,=W, - 2w,+w_, =3(w, +w_,)-2. (3) 


Equations (1), (2), and (3) determine w,, wo, Wx 
and therefore the w, are determined by the speci- 
fication of P, and P,. Consequently the cross sec- 
tion is expressible in terms of P, and P,. From 
(2) and (3) and the inequalities 0 < Wy <1, wy +W4 
<1 it is seen that -1 <P, <+1, the lower and upper 
limits corresponding to (w,, w_,) =(0, 1) and (1,0), 
respectively. Similarly -2<P,<1, the lower lim- 
it corresponding to w,=1 and the upper to w,=0. 

Parallel and perpendicular arrangements of the 
orientation axis relatively to the gamma-ray beam 
direction will be distinguished, respectively, by 
subscripts / for longitudinal, and ¢ for transverse. 
In the longitudinal case the differential cross sec- 
tion for protons with unpolarized y rays is 


0 (6) =a+ bS* + cC +dS*C +eS*C* + P,[-2g - 2hC - (2i + j)S* - (2k + 1)S*C - (2m +n)S*C? - pS*C - qS*C*}, (4) 
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where 
S=sin@, C=cosé. (4’) 


The coefficients are given in Table I. The angle 
between the direction of the ejected proton and 
the gamma-ray beam is denoted by 6. The angu- 
lar distribution is obviously axially symmetric 
in this case and @ suffices for the specification 
of oj. The polarization parameter P, does not 
enter Eq. (4), since the levels .=+1 contribute 
equally to oj, aS may be seen by analyzing the 





incident beam in terms of circularly polarized 
photons. Reflection in a plane through the photon 
beam axis shows in fact that right-handed pho- 
tons incident on 4=+1 contribute equally with left- 
handed ones on p=-1 and vice versa. 

If the orientation axis is perpendicular to the y- 
ray beam, it is convenient to use polar coordinate 
angles (6, #) for the proton direction with the pho- 
ton beam as polar axis and the azimuthal angle $ 
measured from the orientation axis in the right- 
handed sense. The differential cross section for 
protons in this “transverse” case is 


0,(8, &) =a + bS? + cC + dS?C +eS?C? +Pilg+hC + iS? + RS*C +mS*C? + [jS? + 1S?C +nS?C? + pS*C + gS*C?|cos*o} 


+ P,{S[r + sC + tS? +uS?C]sin$}. 


For checking calculations and possibly also ex- 
periments, it is useful to have available the rela- 
tion 


slo, (8, ) + 0 (8, }+7/2)+ 0 (8)] [o(6)] 


0, -2- 0, 0° 


(6) 


The values of P,, P, are here indicated as sub- 
scripts. The left-hand side refers to complete 
alignment and zero polarization, the right-hand 
side to unoriented deuterons. The relation fol- 
lows from the fact that P,=0, P,=-2 corresponds 
to w,=w.,=0, w,=1 which implies that, keeping 


(5) 





values of the incident gamma-ray energy Ey as 
in reference 1. The dependence of 07 (9), o;(@,0), 
oz(6, 7/2) on @ is shown in Figs. 1, 2, 3, respec- 
tively, for E, =22.2 Mev and 102.2 Mev. The 
differential cross section o(@) for the unoriented 
case is shown for comparison by the dashed 
curves when this is feasible. The plots refer 

to three cases I, II, III corresponding to cross 
sections 07, Oy, 7, respectively. Case I cor- 
responds to complete alignment along the orien- 
tation axis and zero polarization, i.e., P,=0, P, 
=1, w,=w_,=4, w,=0. Case II similarly repre- 


the proton direction fixed, the three cross sections sents the condition for the opposite extreme of 


correspond to a complete set of three mutually 
orthogonal spin functions. 

The coefficients a, b, etc., which occur in 
Eqs. (4), (5) are listed in Table I for the same 


the alignment parameter and zero polarization, 

i.e., P,=0, P,=-2, w,=w.,=0, wo=1. The case 
of maximum polarization and alignment along the 
orientation axis is obtained if w,=1, w,=w.,=9, 


Table I. Values of coefficients for Eqs. (4) and (5). 














Ey (Mev) a b c d e g h i j k 
22.2 5.26 51.6 0.844 15.5 1.16 0.83 0.104 0.111 -4,57 0.115 
32.2 §.57 29.1 1.04 10.9 1.00 1.01 0.105 0.329 -5. 76 0.393 
62.2 5.42 7.54 1.05 4.41 0.67 1.45 0.094 0.529 -7.14 0.932 

102.2 4.11 2.50 0.796 2.01 0.401 1.14 0.053 0.426 -5.56 0.82 
152.2 2.87 1.18 0.616 1.02 0.219 0. 766 0.044 0.35 -3.97 0.525 
177.23 2.45 1.13 0.582 0.876 0.175 0.58 0.035 0.411 -3.57 0.455 

Y m n p q Y Ss t u 

22.2 7.10 -0,029 29.9 -8.78 -0.768 -0.906 -3.40 0.812 0. 046 

32.2 6.79 0.021 25.9 -9.30 -0.795 -0.688 -3.29 0.874 0. 065 

62.2 4.63 0.15 16.2 -8.03 -1.07 -0.177 -0.569 0.874 0.178 

102.2 2.98 0.16 9.74 -5.61 -0.885 -0.008 0.405 0. 746 0.222 

152.2 1.92 0.10 6.18 -3. 64 -0.623 0.061 0.561 0.502 0.194 

o2 1.74 0.085 5.80 0. 
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20 40 60 80 40 160 8=— 180 
@ in degrees 


FIG, 1. Differential cross section for the d(y,n)p 
reaction with unpolarized gamma rays of energy 22.2 
and 102.2 Mev and orientation axis being along the gam- 
ma-ray beam. Subscripts I, I, Il correspond to (P;, 
P») having the values (0,1), (0,-2), (1,1), respectively. 


for which P,=P,=1. This condition is referred 
to as Case III. The cross section for any values 
of P, and P, is obtainable from the graphs by 
means of 


o(0, ; P,,P,)=(§+8P, -P,)o,+(4-4P,)o, +P, 0 


1’° 2 2 
(7) 


The differential cross section for unoriented deu- 
terons is obtainable according to this formula as 
o = (20, + 077) /3. The difference between oy and o 
is therefore expected to be larger in absolute 
value than that between oy and o. 

According to (4) and (5), the polarization param- 
eter P, does not affect either o7() or o7(@,0). Its 
effect on o;(0, 7/2) is also small as is seen from 
Fig. 3 and the fact that according to Eq. (7) P, 
enters through (oy - oy) P}- 

The authors are indebted to Professor G. Breit 
for many valuable consultations and particularly 
for having planned the arrangement of the calcu- 
lation’ in terms of amplitudes with a view towards 
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1 1 


O,(8,0) in 





@ in degrees 


FIG. 2. Differential cross section for the d(y,n)p 
reaction with unpolarized gamma rays for orientation 
axis perpendicular to gamma-ray beam and proton dir- 
ection in the plane of these two directions. Other con- 
ventions as in Fig. 1. 
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FIG.3. Differential cross section for the d(y,n)p re- 
action with unpolarized gamma rays for the case of 
gamma-ray beam and proton directions both perpen- 
dicular to the orientation axis. Other conventions 4s 
in Fig. 1. 
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extensions such as that reported on here. 
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K,° DECAYS AND INTERACTIONS* 


D. Luers,t I. S. Mittra, W. J. Willis, and S. S. Yamamoto 


Brookhaven National Laboratory, Upton, New York 
(Received August 29, 1961) 


We have performed an experiment with low- 
energy K,° mesons in the BNL 20-in. hydrogen 
bubble chamber. In this Letter, we wish to re- 
port results of two types: those bearing on the 
structure of the strangeness-changing weak in- 
teraction; and those which give information con- 
cerning K°p and K°p interactions. 

The beam was made by 1250+ 40 Mev/c m~ on 
a CH, target. K,°’s at 47° in the lab were selected, 
and charged particles were removed by a sweep- 
ing magnet before the beam passed into the cham- 
ber, some 110 in. from the target. The direction 
of the K,° beam in the chamber is known to +2° 
but the momentum has a spread from ~ 200 to 
~650 Mev/c, due to associated production of 
both A°K°® and £°K°®, Fermi momentum in the 
target, etc. 

We have made the usual assumptions (consistent 
with our observations) that the only common 
modes of K,° decay into charged particles are’: 


K,°+1*+e* +, 
K,° nt +u* +p, 
K,° +1* +77 +7°, 


which we shall call tev, muv, and m7. 

Under the present conditions, the K,° momen- 
tum is calculable from measurements of the de- 
cay particles but is not overdetermined. (There 
is in general an ambiguity in the momentum 
which is not serious for the present analysis.) 
The most serious backgrounds, mu and ye decays 
in flight, which might be interpreted as neutral 
decays, were eliminated by restrictions on the 
events which also threw out a small and calcula- 
ble fraction of the K,° decays (4 to 6% for the 
different decay modes). 

A sample of 124 identified sev decays was ob- 
tained by identifying all electrons with momen- 


tum <200 Mev/c by bubble density. Any track 
which was not clearly identified on inspection 
was measured by gap counting.” (A muon at the 
limiting momentum has 1/6?=1.29.) The remain- 
der of the events was classified into two groups, 
one containing 52 events which gave a physical 
value of the K,° momentum when interpreted as 

a mmm decay, and a group of 222 events which are 
mév Or TUL. 

If polarizations are not measured, the configu- 
ration of a three-body decay is specified by two 
variables. We choose to use the energy of the 
charged particles in the center-of-mass system 
of the K,°. For each configuration, we have de- 
termined the probability of identifying a tev and 
of not rejecting each of the modes. This involves 
integrating over the uninteresting angles and over 
the experimentally observed K,° spectrum. Given 
this information, we are prepared to make sever- 
al comparisons with the theory. 

1. The ratio 1-et+v/nte~v=1.16+0.17. This 
ratio is one if CP is conserved, but no larger 
than 1.08 in any case.° 

2. Pais and Treiman‘ have shown that the elec- 
tron energy spectrum in the tev decay for a given 
pion energy has a form which is dependent only 
on the type of interaction involved. The pion spec- 
trum is given, to be multiplied by a form factor 
which is expected to be roughly constant.** We 
have divided the Dalitz plot into three strips in 
Tz, each of which is subdivided as shown in 
Fig. 1. 

We multiply the probability for observation 
and detection at each point times the expected 
distribution function for pure scalar, vector, 
and tensor interactions and integrate over each 
region. Theory and experiment are normalized 
to one in regions 1 and 2; 3, 4, and 5; and 6, 7, 
8, and 9; thereby giving a comparison essentially 
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FIG. 1. Regions on Dalitz 

— plot for tev to be used for com- 
parison with experiment are 
shown, together with the distri- 
" bution of events. 





































































































REGION OF DALITZ PLOT 


FIG. 2. Comparison of tev events with theory. 
Histograms represent the theoretical distribution of 
events, corrected for detection efficiency and normal- 
ized in each set of regions with the same range of T;. 
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independent of the form factor. The results are 
shown in Fig. 2. Tensor is excluded; scalar and 


3. If the interaction is assumed to be scalar, 
the form factor is found to increase by a factor 
of ~14 from the region with T, <76 Mev to the 
region with T; >104 Mev. Such a large increase 
with pion energy is not expected, and probably 
allows this interaction to be excluded. 

4. The vector form factor (Fig. 3) is roughly 
constant. The conclusion, then, is that the data 
are consistent with a pure vector interaction-in 
agreement with the universal V-A theory.” 

5. If that part of the form factor due to strong 
interactions is constant, and an intermediate 
boson exists, the energy dependence of the form 
factor may be given in terms of the boson mass, 
as shown by the curves on Fig. 3. From the 
present data a boson with mass >500 Mev/c? 


6. The detection and identification probabilities 
may be integrated over all the decay configura- 
tions with weights appropriate to the vector inter- 
action. This allows the nev branching ratio to be 
determined. The number of tev (0.1%) and mv 
(4.4%) in the sample fitting 777 may be calculated 
This allows a complete set of K,° branching ratios 
into charged modes to be given in Table I. 

7. Given the A/= § rule,® the A/= 4 current 
rule,?° and the K+ branching ratios“ and lifetime, 
the K,° branching ratios may be predicted, Table 
I. There is reasonable agreement with the exper 
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CURVES SHOW THE EFFECT OF 


2.0 AN INTERMEDIATE BOSON _ 

















0 50 100 150 


FIG, 3. The square of the form factor for tev decay 
with vector interaction, The curves show the form- 
factor dependence due to a possible intermediate boson. 


imental numbers. If a AJ = current is assumed, 
together with AS = AQ,*° the numbers obtained, 


also in Table II, are in disagreement with the ex- 


periment. Since there is evidence for A/=4, but 
little information concerning the currents, we 
see this as a test of the current rule. 

8. The 7° energy spectrum in m7 decay has 
been given by Weinberg” and Sawyer and Wali’® 


as W(T 70) =(1+aTqzo)(phase space). Using the AJ 
=} rule, he predicts a= -0.0109+0.022. We find 
a=0.0171+ 0.0065. 

We turn now to the interactions. We have ob- 
served 111 events representing comparable num- 
bers of 


K,°+p ~+K,°+p, 
K,°+p+A°+n", 
K,°+p+D°+m. 


Biswas** has pointed out that a measurement of 
K,°p/(A°n* +22°r*) provides a method for dis- 
tinguishing among the various scattering length 
solutions for Kp interactions, if the K*-N scat- 
tering is known. 

Unfortunately, most of our K,°’s have a momen- 
tum greater than the limit of applicability of the 
zero-range theory (estimated to lie between 200 
and 300 Mev/c).*® The results are given in Fig. 
4. Also shown are shaded regions indicating the 
result which might be expected, within the errors 
on the scattering lengths, using the K scattering 
lengths given by a recent x” analysis of K~-p in- 
teractions.*® Also shown are the range of expect- 
ed values for the ratio at 175 Mev/c, using the 
four Dalitz solutions’” to the same data: The mo- 
mentum dependence is similar to that of solutions 
(1) and (2). Our data would seem to agree better 


Table II. Values of the quantity « = A°/(A°+ =), T=1. 








K momentum (Mev/c) € 
0 0. 570-38° 
~175 0.40+0.03 
~ 230 0.33 +0.14 
~ 360 0.19+0.07 
~ 500 0.32 20.08 





2From K7 -p data, quoted in reference 17. 
From K -p data, reference 16. 


Table I, Experimental branching ratios for charged modes of K,° decays, together with predictions based on the 
Al =4 rule and the A/ =4 and 4] =4 current rules. The ratio 7uv/mev is given because the individual errors are 


correlated, 








Quantity Experiment Al =4 current Al = current 
Tev 48.4 +5 % 40.5 +4.5% 27.4 +3% 
Tv 38.2 +8 % 43.4 +4.5% 29.1 +3% 
mt 13.4 +1.8% 16.1 +1.0% 43.5 +4% 


Tv / tev 0.79 +0.19 














1,06 +0.15 1,06 +0.15 
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FIG, 4. The experimental val- 
ues of K,°/(A°+2=°) compared 
with the approximate values ex- 
pected from K*- and K~-nucleon 
scattering. 
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with solution (1) than with solution (2), or better 
with a~ or b~ than with a* or b+, or, in general, 
with a T =1 scattering length which is not large 
and positive. A more quantitative conclusion 
must await consideration of effective-range ef- 
fects. 

We have also measured the ratio € = A°/(A°+2), 
T=1, using the charge-independent relation 5°r* 
= +n for K°-p, which is a pure T=1 state. The 
results, with some K~-p values,**’"” are given 
in Table II. 

We would like to acknowledge the aid of Dr. E. L. 
Hart and Dr. P. C. Connelly. 
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METHOD FOR DETERMINING THE SPIN OF THE K-17 RESONANCE 


David O. Caldwell* 


CERN, Geneva, Switzerland 
(Received June 12, 1961; revised manuscript received August 31, 1961) 


In order to understand the recently observed! 
resonance in the K-7 system, which is here 
labeled K*, it is essential to determine its spin. 
Recently there have appeared arguments? favor - 
ing a spin of 1 and a proposal’ for determining 
the spin. The appearance of the latter has prompt- 
ed the writing of this note describing a method 
of spin determination which may be simpler ex- 
perimentally and can perhaps provide more infor- 
mation than that of reference 3. It has been tac- 
itly assumed in that paper that the measurements! 
limit the K* spin to 0 or 1 and the isotopic spin to 
1/2, and that because the decay K* +K+7 is seen 
to be rapid, parity is conserved and hence the K* 
has either spin 0 and parity opposite to that of the 
K, or spin 1 and the same parity as the K. Under 
these same restrictions, if one observes the proc- 
ess 

K+He*+K* +He*, (1) 


then the K* must have spin 1, since for the spin 
0 case angular momentum and parity cannot be 
conserved in this reaction. If (1) is not seen, the 
K* spin assignment can be checked by observing 


K+Het+K*+He*+z7, (2) 


and looking at the K*-decay distribution. 

There are two circumstances under which (i) 
will be forbidden and a third reason why it would 
be inhibited, however. Thus if (1) is observed, 
surely (a) the spin of the K * is greater than zero 
(and presumably one) and (b) the isotopic spin of 
the K* is the same as that of the K, which is 1/2. 
If the K* had an isotopic spin of 3/2, which is al- 
lowed by its decay into K+7, then (1) would not 
occur, and hence this point can be checked for a 
nonzero spin K*. 

The inhibition of (1) will be discussed shortly, 
but first let us investigate what could in principle 






be learned from observing the decay angular dis- 
tributions from the K* produced by (1). These 
results will then be applied to (2), which is not in- 
hibited, but which is more difficult to discuss. 

It is well known‘ that one can obtain information 
on the spin of an unstable particle if one observes 
such particles produced near the beam direction. 
A reaction like (1) is, however, particularly use- 
ful for this purpose because the decay angular 
distribution for the K* is uniquely determined by 
its spin, and because most of the observed re- 
actions could be utilized for this analysis. In 
justifying the first of these statements, let us 
begin by considering the initial state of (1): There 
can be no orbital angular momentum (i) about the 
direction (z) of the incident K, and since there is 
no spin, the total angular momentum in that direc- 
tion ( jz) is zero also. Therefore in the final state 
j, =9 too, and since for forward production /, =0, 
then the spin (S) can have no component inz direc- 
tion. Thus for a K* of any possible spin, if the 
K* is produced forward there will be a unique 
angular distribution® for its decay, given by the 
spherical harmonic I¥s o(6)l?, where 6 is the 
angle between the direction of the incoming K and 
the decay momentum in the center-of-mass sys- 
tem of the K*. According to the uncertainty princi- 
ple, the decay distribution should not change ap- 
preciably over K* emission angles of the order of 
the reciprocal of the largest orbital angular mo- 
mentum contributing to the production process, 
which is ~i/pR, where p is the incident K momen- 
tum and R is the “radius” of He*. This, however, 
is just the typical emission angle for events in 
which the He* stays bound. To keep the He* to- 
gether, the momentum transferred to the He*, q, 
must be kept small. From the uncertainty princi- 
ple, g~h/R~0.2 Bev/c. To satisfy this condition, 
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p must be rather large (21.5 Bev/c), and hence 
the emission angle, ~q/p~h/pR is quite forward 
and is just that obtained above as a limiting one 
for the spin determination. Thus most events 
would be useful ones for this spin analysis, which 
not only permits distinguishing between the likely 
values of 0 or 1, but also affords a means to meas- 
ure possibly higher values in a less ambiguous 
way than is presently available.’ 

There is a second method of using (1) [and (2), 
as will be discussed later] to determine the K* 
spin, and this method can be utilized also to check 
parity conservation in the production and decay of 
the K*, One uses not-too-forward events so as to 
define better the direction, x, of the normal to the 
plane formed by the relative momentum of the in- 
cident particles and that of the outgoing particles 
(before K* decay). Performing an inversion about 
that plane followed by a 180° rotation about x would 
change the initial wave function by (-1)“ K, where 
Px is the parity of fhe K, and the final wave func- 
tion by (-1)?K (-1)"*.° Therefore, if Py*=Py, 
the projection of the k* spin in the x direction, 
S,, is even, whereas if Px*=-Px, Sx is odd. 
Thus for example if the K* spin is 1, the decay 
distribution should be given by |Y, ,(6’)I? < cos*@’ 
(where 6’ is the angle between x and the decay 
momentum in the rest system of the K*), because 
Px*=Px for S=1 if parity is conserved in 
K*+K+m,. However, if parity is not conserved, 
then Py* can be -Px and one would see a sin*6’ 
distribution. 

Thus one can in principle determine from ob- 
servations of (1) the spin of the K*, no matter 
how large it is, and also get information on its 
isotopic spin and on whether parity is conserved 
in its decay or production. However, the previ- 
ously mentioned inhibition of this reaction limits 
its usefulness in practice. That (1) will not occur 
with large probability even if it is allowed by spin 
and isotopic spin considerations can be seen by 
considering the case in which the K * is produced 
forward, followed by K* decay also along the beam 
direction. In this case one has K + He*+K+7+ He*, 
with all the momenta essentially along a line, the 
z axis. If parity is conserved in the production 
and decay process, a reflection of the system 
about any plane containing the z axis should not 
alter the parity. However, this is not true for the 
postulated system because the 7 is pseudoscalar, 
and there can be no compensating parity change 
provided by angular momentum, since there is no 
orbital angular momentum about the z axis. Since 
we have seen above that forward K* decays would 
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otherwise be highly probable, this conflict can be 
resolved only if the amplitude for K : production 

goes to zero in the forward direction. Thus re- 

action (1) is inhibited, and it would be difficult to 
obtain the information for the analysis of the de- 
cay distributions. 

This analysis can, however, be applied to re- 
action (2), which is not forbidden by spin or iso- 
topic spin considerations, and is not inhibited as 
(1) is, because for (2) there are two pions in the 
final state, and hence there is no difficulty in con- 
serving the parity. For (2), one must look at those 
cases in which both the K* and the 7 travel near 
the z direction; then, as with (1), the K ” decay 
distribution about z is the same unique function of 
the K* spin. Even the analysis of the decay distri- 
butions with respect to x to determine spin and to 
check parity conservation can be applied to (2) 
provided x is as defined above; i.e., the K*, 7, 
and He* lie essentially in a plane, although the 
noncoplanarity angle can again be of the order of 
h/pR. Then because of the additional 7, Py*=Px 
makes S, odd, whereas Py*=-Px makes S, even, 

An important experimental problem in any of the 
above observations is being sure that one has He 
in the final state. Fortunately He* has no excited 
states, and its disintegration into hydrogen iso- 
topes is clearly distinguishable, but the emission 
of He*+n instead of He* provides some difficulty. 
Of presently available techniques, the helium 
bubble chamber appears the most suitable for 
this investigation, and in such a device with a 
magnetic field the He*+n events should usually 
be rejected on the basis of a kinematic fit. In 
some cases, the fraction of the total depending 
strongly upon the incident K energy and the 
chamber size, the He*’s and He*’s can be separa- 
ted on the basis of momentum and range. In fact, 
observation of the He* range and angle alone clear- 
ly separates reaction (1) from other processes, 
such as elastic scattering. Reaction (2) is harder 
to measure than (1), but all of the final particles 
can be charged, and hence clear cases can be ob- 
tained. Thus although the mere observance of (1) 
can determine whether the K* spin is 0 or 1 under 
the highly probable assumptions of reference 3, 
it seems possible to settle this important question 
with fewer restrictions by observing the decay 
distributions from (2) and perhaps from the much 
rarer (1). 

I wish to express my gratitude to B. T. Feld 
and D. B. Lichtenberg for discussions, and par- 
ticularly to T. D. Lee, who considerably altered 
the contents of this note. 
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We have studied the reactions Li®(p, 2p)He® and 
Li"(p, 2p)He® by detecting the product protons in 
coincidence. The 155-Mev proton beam of the 
Orsay synchrocyclotron had a 10-20% duty cycle.” 
The experimental apparatus has been used for a 
similar study in carbon.’ It consists of two tele- 
scopes, each with (1) two plastic scintillation 
counters to define the scattered beams and (2) 
one NaI counter to measure the energy of each 
product proton. The two telescopes are used to 
study the angular distribution of two protons emit- 
ted in coincidence in a plane which includes the in- 
cident beam and at angles @, and @, with the inci- 
dent beam and arranged symmetrically around it. 
The angular resolution in the scattering plane is 
+2.5°. At each angle, the summed energy of the 
two coincident protons is displayed on a pulse- 
height analyzer when the energy of one of the two 
protons lies between 62 and 82 Mev for Li® and 60 
and 80 Mev for Li’. This limits the events studied 
to an energy ratio, x?=E,/E,, of 0.7<x?<1.0 for 
low-binding-energy protons and 0.6<x?<1.0 for 
high-binding-energy protons. The energy cali- 
bration is carried out with free p-p scattering 
events from a CH, target. The targets were pre- 
pared from material supplied by Oak Ridge Na- 
tional Laboratory and had a purity of 95.6% (Li®) 
and 99.9% (Li’). 

The summed energy spectra at 6, =6,=40° dis- 
played in Fig. 1 are typical. Two peaks are ob- 
Served for each target, corresponding to binding 
energies of the target proton of 23.54 0.7 Mev and 
10.0: 1.4 Mev for Li’ and 20.0+ 0.7 Mev and 4+1.4 
Mev for Li®.* The position and width of the low- 
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binding-energy peaks correspond to the residual 
nuclei He® and He® being left in their ground state, 
although excitation of the 1.71-Mev level* in He® 
cannot be excluded. The binding energies involved 
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FIG. 1. Summed energy spectra of Li® and Li’ at 
6,=0,=40°. The abscissa also indicates the binding 
energy of the target protons. 
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FIG. 2. Angular correlation 
distributions for Li’. 
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suggest that the higher binding energy peaks are 
produced by the interaction of the incident protons 
with “inner core” protons in the target nuclei 
while the lower binding energy peaks are pro- 
duced by the interaction with “outer shell” pro- 
tons. This conclusion is based on the fact that 
the interaction time in the (p, 2p) reaction is 
much smaller than the rearrangement time of 
the residual nucleus. Even for the ground state 
of He’, these two times differ by a factor of 30 
as calculated from the level width of 0.5 Mev re- 
ported in reference 4. The binding energies for 
the inner core protons are similar to those for 
He*. 

We have obtained angular correlation distribu- 
tions from energy spectra measured at several 
angles 6 =06, =6, by integrating the number of 
events observed in a fixed summed energy inter- 
val of 15 Mev centered around each peak. For 
the higher binding energy protons this method 
yields only an upper limit to the cross section 
at small and large angles (@<30°, @>50°) since 
the corresponding peaks ride on a continuum 
which becomes appreciable. The results are 
presented in Figs. 2 and3. The cross-section 
scale is in arbitrary units, but the same units 
are used for all curves. The errors shown are 
of statistical origin. The results for Li’ are in 
agreement with previous studies® with ordinary 
targets in the regions where comparison is pos- 
sible. 
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I 
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The angular correlation distributions for Li’ 
(Fig. 2) have characteristic features previously 
observed with C**.?"5** The high-binding-energy 
peak has a maximum at about 41° as expected 
for s-shell protons, and the low-binding-energy 
peak has a minimum at about 42°, characteristic 
of p-shell protons. At small and large angles 
the cross sections fall off much more rapidly than 
for C’*, On the usual Born approximation inter- 
pretation, this observation means that the wave 
functions of the target protons in Li’ have fewer 
high-momentum components than in C™*, The 
same conclusion has been reached from a study 
of the energy-sharing’ in the same reaction. 

The angular correlation distribution for the 
high-binding-energy protons in Li® (Fig. 3) is 
very similar to that of the corresponding pro- 
tons in Li’; the 23% difference in peak cross 
section can be explained by the difference in ab- 
sorption. However, the angular correlation dis- 
tributions of the low-binding-energy protons in 
Li® and Li’ are very different, that for Li® being 
instead very similar to that of the high-binding- 
energy protons. The observation of the narrow 
peaks in the summed energy spectrum suggests 
that this difference is not the result of the differ 
ence in lifetime of the residual nuclei but rather 
due to differences in the properties of the low- 
binding-energy proton in Li® and Li’. Using the 
impulse approximation interpretation, the wave 
function of the low-binding-energy proton in Li’ 
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contains many more low-momentum components 
than the corresponding proton in Li’. In fact our 
Li’ distribution is consistent with an s-proton 
interpretation, but we do not know whether a p 
proton with important low-momentum components 
might not yield a similar distribution. 

The observed dip at about 41° in the angular 
correlation distribution for the low-binding-ener - 
gy protons in Li’ indicates that these protons have 
a pronounced minimum in their linear momentum 
distribution at zero momentum. Therefore the 
shell-model description of the low-binding-energy 
protons in Li’ as p-shell protons is consistent 
with our observations. The results with the low- 
binding-energy protons in Li® are difficult to un- 
derstand on the basis of the usual shell model, 
especially the large difference between the Li® 
and Li’ distributions. Target protons with a max- 
imum in their linear momentum distribution at 
zero momentum yield an angular correlation dis- 
tribution with a maximum at about 41°. A model 
for Li® which yields such low-binding-energy s 
protons (a description consistent with our results) 
is the (alpha + deuteron) complex, in which the Li® 
spin is entirely due to the intrinsic spins of the 
neutron and proton in the deuteron. 

The differences reported in this note seem to 
Ws to illustrate well the usefulness of studies of 


30 





40 60 6 


(p, 2p) reactions at high energies in order to ob- 
tain information on nuclear structure. A com- 
plete report of similar studies with light ele- 
ments will be published elsewhere. 

We wish to thank the synchrocyclotron crew 
and the electronics group for their assistance. 
One of us (K. S.) wishes to thank the Laboratoire 
Joliot-Curie de Physique Nucléaire for their hos- 
pitality, and the Alfred P. Sloan Foundation for 
support. 
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The purpose of this note is to report a result for 
the helicity of the proton from A decay. We find 
that our data strongly favor positive helicity. 

This result is obtained from a study of about 
2000 A’s produced by K~ absorption in He using 
the Duke helium chamber,' and a low-energy K~ 
beam at Berkeley.? The study, based on 105 pro- 
ton scattering events from A decays, was made 
as part of a continuing experiment, a portion of 
which is described in detail elsewhere.*® 

The helicity of the proton in A decay is defined 
as the sign of the proton polarization measured 
in the decay rest frame. As shown by several 
authors,* this polarization is given by 


3 = “a k*/ik"|, 


(1) 


where a is the asymmetry parameter of the decay 
A+m1~+p, or in terms of the s and p amplitudes, 
a = 2Re(s*p)/(s? +p?) and k* is the c.m. momentum 
of the proton. The proton polarization, longitudi- 
nal in the A rest frame, becomes partially trans- 
verse when transformed into the laboratory frame 
since (nonrelativistically) the spin remains rigid 
while the momentum k* transforms into k,(lab), 
making an angle © with respect to k*. The trans- 
verse component, osinl, can then be measured 
by means of subsequent elastic scattering. 

In principle both the magnitude and the sign of 
a@ can be obtained by measuring the (azimuthal) 
angular distribution of the scattered protons. In 
practice, however, since only a small portion of 
the available data has as yet been analyzed, only 
the sign of a is determined with good accuracy in 
this report. 

Severe limits on the magnitude of a have already 
been set from the results of the up-down asym- 
metry experiments at Brookhaven and Berkeley.® 
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The latest unpublished Berkeley results® indicate 
that |a| 20.85. We employ this result as a con- 
straint on the determination of the sign of a in 
the present experiment. Details of the analysis 
are given below. 

If an incoming proton beam of (lab) momentum 
k; and arbitrary polarization @ is scattered from 
a spin 0 nucleus, the angular distribution’ of the 
elastically scattered protons which emerge with 
polarization P,,,(6,€) is 


f(0,€)~1+5-P_ (6, €), (2) 





where @= (proton) center-of-mass scattering angle, 
€=lab proton energy at the point of scattering, 
k= outgoing proton lab momentum, and Pan(,€) 
=P, (6, €)k; xK,/ Ik, xKy| = analyzing polarization’ 
This expression can be written in a form more 
useful for analysis by expressing o Pan(6, €) in 
terms of the angle [ defined above, and ¥, an 
azimuthal scattering angle which is conveniently 
chosen as the angle between the A decay plane 
and the p=-He scattering plane. The distribution 
then becomes 


f(¥)dQ=[1+a sin P (6, €) sin’ ]dQ. (3) 





The quantities [ and W are directly obtainable’ 
from a kinematic analysis of each event, and sinc? 
the analyzing polarization is known?® (see Fig. 1), 
a may be determined by an appropriate best-fit 
procedure using (3) as the theoretical distribution. 

The important characteristics of the various 
subsamples chosen for analysis are given in Table 
I. Of the total of ~ 150 A-decay proton scattering 
events found by the scanners, only 105 survived 
the selection criteria for sample A which were 
imposed in order to remove the Coulomb scatter 
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FIG, 1. Analyzing polarization for proton helium scat- 
tering as a function of c.m. angle for various proton (lab) 
kinetic energies. Data at and above 10 Mev are taken 
from J. L. Gammel and R, M, Thaler, Phys. Rev. 109, 
2041 (1951). Data below 10 Mev are calculated from 
phase shifts given by A, Juveland and W. Jentschke, 

Z. Physik 144, 521 (1956). 


ings and possible inelastic contamination. These 
selection criteria are: (1) @=26°, (2) 4 Mev<e 
<22.5 Mev. Of the total of 105 events in A, 100 
stop in the chamber. The energy spectrum and 
angular distribution of all events is shown in Fig. 
2. 

Each event in the sample was subjected to a full 
kinematic analysis which yielded best-fit values 


C.M. SCATTERING ANGLE (DEGREES) 


FIG, 2. (a) Proton energy spectrum of Sample A 
(02 26°, 4 Mev<€<22.5 Mev). (b) Center-of-mass 
angular distribution of Sample A. 


of the dip and azimuth of each track. This infor- 
mation, when coupled with the accurate energy 
determination provided by the range of the scat- 
tered proton, permitted the calculation of sin¥ 
and P,, to an average accuracy of 15% and 10%, 
respectively." 

The determination of the sign of a was carried 
out using three techniques, (a) the method of maxi- 
mum likelihood, (b) comparison of the measured 


Table I, Characteristics of the proton scattering sample. N=No. of events in sample; €=average proton energy 
at scattering point; @= average c.m. scattering angle; Pg,=average analyzing polarization; (P,y Sin!) ay = average 
effective polarization of sample; 7. experimental value of right-left asymmetry; (n)= predicted value of right- 

P 


left asymmetry using the measured (Pan sinl)ay. 





Sample definition N € 


) 8 


Pp 4sinT) (P_sint) 7 
an av an av exp 





A. 0>26°,4 Mev<e<22.5 Mev 
B. 6>26°,4 Mev<€<22.5 Mev, P Sint > 0.25 
C, @>40°, 4 Mev<e<22.5 Mev 


=... 


105 12.0 0.36 
57 10.6 0.47 
85 11.8 0.41 


0.77 0.28 0.20 -0.18a@ 60.5° 
0.86 0.40 0.15 -0.26a 69.0° 
0.76 0.31 0.20 -0.20a 66.5° 
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right-left asymmetry with its expected value, and 


(c) a x? fit of (7) to the © distribution. 

The likelihood function, L(a), representing the 
total probability for observing our 105 events in 
sample A, is a product of the distributions (3) 
evaluated for each event: 


L(a)= piae - oP, 6;, €) sinl sin¥ (4) 
A plot of L(a) vs a for sample A is shown in 


Fig. 3. Although the best-fit magnitude of a, as 
determined from the position of the maximum of 


L(a), is quite consistent with the limits set by the 


up-down asymmetry experiments (0.85 < |a|<1), 
it is clear from the width of L(a) that the present 
data cannot determine |a@| with comparable pre- 
cision. We therefore accept the previous deter- 
mination of |a|~1, and ask the question: “Given 
the magnitude of a, what is the likelihood of the 
hypothesis a=~-1, relative to the hypothesis 
a@=+1”? One such relative likelihood measure 
is given by the ratio L(a= -1)/L(a=+1)=2x10° 
which strongly favors” a= -1, 

Another such relative measure can be obtained 
by comparison of the measured azimuthal angle 
asymmetry 7, defined’* by 


" [NV , (sin¥ 20)-N  sinw <0) VW, +N) (5) 


with its expected value (n) given by 
(n) = (2/m)a(sinl P,(0,€)),.- (6) 


For our sample A, (n)=-0.18a, and the meas- 
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FIG. 3. Likelihood function L(@) vs a for Samples A, 


B, andC. The shaded areas represent regions allowed 
Sample B is defined 


by up-down asymmetry experiment. 


ured asymmetry is Noxp=9.20+0.1. Thus the 
data are consistent with the requirement |a|~=1, 
and are in good agreement with the hypothesis 
a=-l, but are in disagreement with the hypothe- 
sis a=+1 by about 3.8 standard deviations,” 

An alternate approach is based on a x? analysis 
of the observed distribution in © (see Fig. 4), 
The expected distributions g(¥)d¥ for a=1, re- 
spectively, follow from (3), integrated over solid 
angle: 


g(¥)d¥ = [1+ asin’ P, €)) ay sinv|]dv, (7a) 
g()d¥ = [1 -a(sin? P, €)) ay sin’ |d¥, (7b) 


where (sin P,,,(0,€))gy is the experimental aver- 
age given in Table I.* The most probable value 
of x? for this fit is 5.0. The distribution (7a) is 
in excellent agreement with the data of Fig. 4, 
yielding a x? of 3.9, corresponding to a probabili- 
ty of 56% for the observation of the data, if a=-l. 
On the other hand, the distribution (7b) has a y’ 
of 25.7, corresponding to a probability of less 
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FIG, 4. Distribution of angle ¥ between A decay and 


by P,, sint > 0.25, 4 Mev<e€<22.5 Mev; sample C is de- proton scattering planes. The expected distributions for 


fined by @> 40°, 4 Mev<e < 22.5 Mev. 
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@ = +1 are represented by smooth curves. 
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than 1% for the observation of the data, if a=+1. 

Aside from the purely statistical errors given 
above, we have considered several possible 
sources of systematic error—inelastic contamin- 
ation, Coulomb contamination, scanning bias 
against events in certain regions of energy and 
scattering angle (see Fig. 2), and finally polari- 
zation of the A sample. We believe all of these 
effects to be either small or nonexistent for rea- 
sons discussed below. 

Firstly, each of the scattered protons had an 
energy below threshold for p+ He*~d+He*, which 
constitutes 90% of the inelastic cross section at 
low energy. Secondly, the @ = 26% criterion was 
designed to eliminate ~ 90% of the Coulomb events. 
In any event a small contamination of Coulomb 
scatters will not affect the maximum-likelihood 
estimate of a because these events, at small 6, 
would be assigned a small Pan (see Fig. 1). 

In principle, since the angle v between the 
production and decay planes is uncorrelated with 
scattering energy and the angle, the bias against 
low € and small 6 events should not affect the v 
distribution. To test this point experimentally, 
we selected two subsamples B and C according 
to the criteria given in Table I. B is a high ef- 
fective polarization sample rich in low-energy 
events; C is a sample composed of large scatter- 
ing angle events. Both the likelihood functions 
(see Fig. 3) and the asymmetries (see Table I) 
for B and C are in good agreement with their ex- 
pected values, thus confirming the internal con- 
sistency of the data. 

Lastly, we consider the possibility of the exis- 
tence of a polarization of the A sample, which 
would complicate the analysis to some extent.* 
The majority of the A’s in the sample were pro- 
duced via a two-step © -conversion process,’® 
ie., 


K~+N-Z+17, (8a) 
D+N’?A+N’. (8b) 


The outgoing A’s, if polarized at all, must be 
polarized perpendicular to the production plane 
of (8b). However, since all such production 
Planes are accepted, the outgoing A’s must be 
polarized, 

Owing to all the above considerations, we in- 
clude no systematic effects in our estimation of 
the errors in this experiment. 

In order to make a direct comparison of our 
Tesults with those of other experiments of this 
type, we must disregard the constraint placed on 
the magnitude of a by the up-down asymmetry ex- 


periments. Then, from the likelihood function of 
Fig. 3, the most likely value of a is 


a = 0.75 70035, (9) 


where the error on the positive side corresponds 
to the half-width of L(a). The result (9) can be 
compared with those of Boldt et al.*® and Birge 
et al.,’” who find a =+0.85+ 0.20 and -0.45+ 0.5, 
respectively. Our result is clearly in disagree- 
ment with the former; it is in agreement with 
that of Birge et al., but is even more indicative 
of a negative a. 

Furthermore, we should like to emphasize that 
by taking the known limits on |a| properly into 
account, as described above, one obtains consid- 
erably stronger results. We prefer to summarize 
these results as follows: For values of |a| con- 
sistent with the known limits, our data favor the 
hypothesis that a is negative over the hypothesis 
that a is positive by at least 3 standard deviations. 

One of us (JL) would like to acknowledge some 
stimulating correspondence with Dr. R. Birge. 
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ERRATUM 





SELF-CONSISTENT CALCULATION OF THE 
MASS AND WIDTH OF THE J=1, T=1, 77 RESO- 
NANCE. Fredrik Zachariasen [Phys. Rev. Let- 


ters 7, 112 (1961)]. 


Due to an error in the numerical computations 
the mass value mp ~ 950 Mev is wrong and should 
be changed to Mo ~350 Mev. While the agreement 
with the experimental value is slightly less good, 
it is perhaps encouraging to have too large a 
coupling constant going with too small a mass. 

An improved calculation might then give a small- 
er coupling constant, corresponding to a weaker 
attraction which could be consistent with a larger 
mass, while it is hard to see how a smaller coup- 
ling constant could go with a smaller mass. 

I am indebted to J. Mathews for pointing out 
that the mass was incorrect, and to C. Zemach 
for aid in computing the correct value. 
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